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THE INTERACTION OF MODERATELY STRONG SHOCK WAVES 
WITH THICK PERFORATED WALLS OF LOW POROSITY 

by 

Daniel J. Grant 
Goddard Space Flight Center 


INTRODUCTION 

The present work is concerned with a means of theoretically predicting the interaction of a travel- 
ing plane shock wave with a thick perforated wall. The flows produced within the perforations, after 
the plane wave reflects from the upstream face of the plate, are influenced by the combined effects of 
heat transfer to the wall and friction with the wall. These effects attenuate the strength of the trans- 
mitted wave and, therefore, enable the plate to function as a wave suppressor. These same effects 
were first considered in the quenching of detonation waves in coal-gas/air mixtures. 

Systematic research of the phenomena of flashback (unstable combustion in which the flame 
front propagates back through a fuel feed device and into the fuel supply system) in combustible gas 
mixtures has been reported as early as 1 8 1 5 by Sir Humphry Davy (Reference 1 ) in his investigation 
of coal mine gas explosions. His experiments on the quenching of coal-gas/air combustion waves by 
small-diameter tubes and fine mesh wire gauze resulted in his invention of the miner’s safety lamp. He 
was one of the first to recognize the importance of the ratio L/d and heat transfer in suppressing flash- 
back. Over one century later, Pay man and Wheeler (Reference 2) investigated the effects of flame 
speed, tube conductivity, and mixture ratio on quenching distances by means of similar wave- 
attenuation experiments. Later work (References 3 to 5) further emphasized heat transfer and throt- 
tling aspects of narrow channels intended as flashback arrestors in burners or torch tips. 

Flashback-type instabilities can occur in rocket chambers using premixed gaseous propellants. 

The thermodynamics and fluid dynamics of these devices are further complicated by the chemical 
reaction occurring behind the detonation waves plus the complex perturbations of the flame velocity 
resulting from interactions with pressure waves reflected from chamber walls and injector faces (Refer- 
ence 6). 

The simplified theoretical model used in the present work considers the flow through an idealized 
injector plate (a metering device that introduces and controls the flow of fluid propellants to a reac- 
tion chamber) produced by the impingement of a traveling shock wave in a nonreactive gas. Such a 
model is necessary in order to predict the attenuation factors for pressure and temperature resulting 
from interaction of the incident wave with the face of the injector plate. Combining the results from 
this study with combustion kinetics and considering the transient heat balance, a numerical 
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computational approach may be developed for the design of injectors for pulse-mode reactors using 
premixed gaseous propellants such that a specified duty cycle can be used without experiencing flash- 
back through the injector to the supply line. 

The attenuation of traveling waves by structural elements has been examined by previous re- 
searchers using shock-tube facilities. The difficulty of obtaining a rigorous theoretical approach to the 
problem can be appreciated when the complexities and ambiguities associated with shock impinge- 
ments on a simple area convergence are considered (Reference 7). As a consequence, most studies 
have been experimental, emphasizing the facility design and instrumentation techniques and primarily 
seeking fundamental theory on wave dynamics in selected flow channels. 

Other investigators have used simplified analytical approaches with varying degrees of experimen- 
tal correlation. Glass and Patterson (Reference 8) demonstrated the slowing down of reflected shock 
waves by screens using a schlieren streak record. Bowman and Niblett (Reference 9) determined wave 
velocities from schlieren photographs of shock-tube plane waves acting upon a woven wire grid of high 
porosity for a single shock strength ( S = 1.65 in air). Dosanjh (Reference 10) used a hot-wire-triggered 
chronograph to measure incident-, reflected-, and transmitted-wave speeds. Shadowgraphs of the 
reflected and transmitted waves were taken for grids having 75 percent porosity ratio. The primary 
purpose of the work was to study shock-wave interaction with a turbulent-flow field produced by a 
grid in a shock tube. The shadowgraphs provided definitive evidence of auxiliary transmitted-wave 
formation resulting from base-region effects accompanying choked flow through the grids. 

Barthel (Reference 11) was concerned with the attenuation of blast waves by rigid grids to pro- 
tect structures. He proposed a simplified mathematical flow model to consider theoretically the type 
of flows that Dosanjh had studied on an empirical basis because of the complexity of the field dis- 
charged from the grid. Barthel’s mathematical model was based on the consideration of the overall 
wave pattern and the local flow pattern in the vicinity of his constrictions. , At some distance from his 
grid, the flow was assumed to be represented by one-dimensional wave-interaction concepts. Flow 
through the grid was based on compressible, one-dimensional steady flow. By considering the com- 
patibility conditions of flow fields behind and in front of each plane wave generated as a result of 
shock impingement with the model, a trial and error method was used to compute the strength of the 
reflected and transmitted waves. For Barthel’s analysis, one-dimensional isentropic flow was assumed 
to exist in the grid. Since his grids were very thin, the residence times of the grid flows were low and 
the assumption of adiabatic conditions was not too unreasonable. In the Mach number range 1.0 to 
2.2, he obtained correlation between predicted and experimental reflected-wave overpressure ratios 
for round-hole grids with sharp edges having 40 and 72.3 percent open area; however, on prediction 
of the transmitted overpressures, the correlation was poor in the Mach number range of 1.6 to 2.2, or 
where choked conditions existed in the grid. This work of Barthel, Dosanjh, and the others (Refer- 
ences 7 to 1 1) was performed with passage length-to-diameter ratios (L/d) of 1 or less. 

Crist (Reference 12), concerned with blast wave attenuation in ventilator systems using rock-bed 
filters, showed that attenuation ratios became independent of bed length when the length was in- 
creased to twice the bed diameter. A lesser dependence on void volume ratio was also demonstrated 
in limited experiments with simple area convergence of the ventilator shaft, using L/d ratios of 8 and 
24 and a fixed length of 4 ft; improved attenuation was evidenced with the smaller openings. This 


i 


2 


I 


situation is closely allied to the case of an injector plate with orifice L/d values of about 10: 1 and 
open area ratios of 25 percent or less. 

The present study makes use of pertinent experimental techniques and theoretical considerations 
to formulate a mathematical flow model for moderately high wave speeds and thick walls with low 
porosity. 

THEORETICAL BACKGROUND AND FLOW MODEL 

It is assumed that the conditions produced by the reflection of the incident shock wave yield 
upstream reservoir conditions that will supply the flow through the idealized injector for a period of 
time sufficiently long to approximate steady-flow conditions. Rudinger (Reference 1 3) and others 
have shown that shock tubes produce appreciable periods of steady flow through short-duct elements 
for low incident-wave Mach numbers. With strong shock waves, the transient periods required for gas 
transport through the longer injector models and for base-pressure adjustment are still only a small 
fraction of the total steady period that will exist before the reflected incident wave collides with its 
upstream contact surface or reflected expansion wave. 

The effects of heat transfer to the injector plate and turbulent friction in the injector plate will 
be considered in the calculations to predict the reflected and transmitted waves and the conditions of 
state behind them. In the range of incident-wave speeds selected, the pressure ratios across the model 
faces will be such that choked flow will exist in every case. A quasi-steady one-dimensional approach 
will be used to define regions upstream from the perforated wall and to define flow through the con- 
stricted passages and in the downstream regions. Base-region flows will also be evaluated in terms of 
steady-state conditions and will be the final relatively steady conditions that will exist after the brief 
interval of time required for secondary flows in the base region to stabilize. The flow in the base 
region can either lower or raise the base-region pressure with respect to initial conditions depending 
upon incident-wave strength; this transient phenomenon will be referred to as base-region pumping. 
This study is confined to incident-wave Mach numbers of 2 to 5, open area or porosity ratios of 25 per- 
cent, and L/d of approximately 10. Dry nitrogen gas will be used as the fluid medium. 

As has been observed with screens and thin gridlike structures (References 9 to 11), the complex 
wavelet pattern, which is established immediately after wave impingement, coalesces at some distance 
from the grid to form a one-dimensional system. With thick perforated walls of high solidity and 
moderately high incident-wave Mach numbers, similar wave patterns are expected with increased tur- 
bulence in the contact flows and adjustments in base flow as a consequence of the higher stagnation 
pressures and greater base areas. 

Therefore, the basis for examination of the specific wave patterns for these flows can be very 
similar to the approach suggested by Barthel (Reference 1 1); however, the flow through the thick 
grids cannot be considered isentropic. The higher range of incident-wave strengths will produce higher 
stagnation temperatures in the reservoir region. The longer path through the injector model will en- 
courage heat transfer from the shock-excited gases. The combination of high thermal gradients, con- 
vective cooling, and entropy gains caused by friction during flow through the model requires that a 
nonisentropic approach be considered. Since the relatively steady pressure ratios across the model 
will be lower than critical, choking is assumed to exist in all cases. This analysis can also be 
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accomplished on a one-dimensional quasi-steady basis by insuring compatibility of conditions behind 
the reflected incident wave and entry conditions at the upstream face of the model. It can be similarly 
accomplished by considering compatibility of final quasi-steady conditions on the downstream face of 
the model and the state of the gas upstream of the transmitted waves. 

Wave Patterns 

Since the incident-wave Mach number range is high enough to cause choking at the entrance or 
exit of the thick-wall grid in all cases, final supersonic regions will appear downstream from the throat 
of the choked region as a consequence of the sudden downstream expansion. Depending upon the 
incident-wave strengths, auxiliary left-facing wave systems will be generated in addition to the trans- 
mitted wave and contact surface. The conditions upstream of these auxiliary waves must be compati- 
ble with the base flows emanating from the perforated plate; the conditions downstream of the 
auxiliary wave must be compatible with conditions downstream of the transmitted-wave contact 
surface. 

Figure 1 indicates two possible wave patterns expected with a range of incident shock strengths 
of 4.5 to 29 (i.e., M = 2 to 5). With a primary ratio equal to the breakoff value (or the value at which 
the streams exhausting from the perforations completely fill the channel), a standing, left-facing 
auxiliary wave should appear (Reference 14). With additional pressure adjustment in the base region 
of the plate, the auxiliary transmitted wave will be modified to become a left-facing traveling shock 
wave (Figure la) as reported by Dosanjh (Reference 10) and Barthel (Reference 1 1). Barthel also 
speculated on the likelihood of a left-traveling rarefaction wave (Figure lb) with further increase of 
incident-wave strength. Although this represents a theoretically compatible condition, this phenomenon 
was not observed in the experiments he conducted. At the higher range of shock strengths, this should 
be considered. 

Contrary to Barthel’s case, with the area ratios under consideration (e = model free area/channel 
area = 0.25) and the incident-wave Mach number range of 2 to 5, the reflected shock wave will travel 



(a) (b) 

Figure 1-Wave patterns for choked flow through perforated walls, (a) Auxiliary traveling shock, (b) Auxiliary 

rarefaction wave. 
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upstream against the drift flow of the incident wave and should increase in velocity as the incident- 
wave strength is increased. Thus, the reflected wave will not be swept downstream to be stabilized by 
the perforated plate as occurred in Barthel’s work. 

The incident wave is assumed to have developed into a plane wave of constant strength by the 
time it reaches the model since the length of the driven section of the shock tube is several times the 
minimum specified 10 widths found to be necessary (Reference 15). 

Influence Coefficients Method 

The incident shock wave, when reflected from the upstream face of the model, will produce a 
reservoir region of high-pressure, high-temperature gas upstream of the model. The Mach number in 
this zone will be subsonic with flow directed toward the model. The gas will enter the sharp-edged, 
cylindrical holes in the model and be subjected simultaneously to cooling and friction. 

It is reasonable to assume that a turbulent boundary layer will develop on the walls of the 
perforations since the reservoir gas has been shocked by both the incident and reflected waves and 
because the sharp-edged entries of the model perforations act as boundary-layer trips. The transient 
mixing of shocks and expansions at the upstream face of the model prior to the establishment of 
steady flow through the model is another contributor to high initial turbulence levels at the entries to 
model ducts. Analysis of the temperature increase of the model wall during the brief period of 
relatively steady flow indicated a negligible temperature rise in the perforated plate. This negligible 
increase and the fact that the plates are made of materials having good thermal conductivity provide a 
reasonable basis for considering the wall temperatures constant during the testing time. 

In the selected range of incident-wave Mach numbers, the passages will choke and therefore deter- 
mine the entrance Mach number. This entrance velocity will in turn establish the particle velocity 
downstream of the reflected wave so that compatible quasi-steady conditions can exist in reservoir 
region 3 of Figure 1 . 

Shapiro (Reference 16) indicated the general analysis and methods of numerical integration 
required for the solution of the working differential equations in terms of influence coefficients for 
generalized one-dimensional steady flow having constant area, constant specific heat, constant molec- 
ular weight, constant wall temperature, and combined friction and heat transfer. To simplify the 
differential equation for heat transfer expressed by means of the energy equation in terms of stagna- 
tion enthalpy changes, the following additional assumptions were made: (1) The adiabatic wall 
temperature T does not differ very much from the stagnation temperature T 0 ; i.e., the recovery 
factor r is unity; and (2) Reynolds’ analogy between friction and heat transfer is valid. These assump- 
tions permit the energy equation to be written as 

dT Q dx 

= 2 /— , ( 1 ) 

T W ~T, d 

and the influence coefficients equation becomes 
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dx 


dM 1 — Ft q 


dT n 


+ FAf — 
d 



( 2 ) 


Since T w is independent of x for a constant wall temperature and considering any pair of sections 
1 and 2, integrating Equation 1 yields 


1 4 /( x 2 - x ,) T w -T 0i 

- = In- 

2d T w - T 02 

T J T oi - 1 

= In 

T w /T 0l ~ T Q2 /T 0l 


(3) 


Thus for any initial value of T 01 and for selected values of f, d, and T w , the variation of the stagnation 
temperature with displacement through the tube or duct can be computed. Equation 2 can then be 
approximately integrated in a finite difference form over the short interval between incremental 
sections to provide the variation of Mach number with distance in the duct. The coefficient of dT Q is 
assumed constant at its mean value during the interval. Fr 0 and F^ denote the values of the influence 
coefficients when evaluated at M = (M l + M 2 )/2 and T Q = ( T 0l + T Q2 )/ 2 for the interval. Then the 
approximate integration for the case of constant wall temperature yields 


M\ - M ] 



Ft o 


+ 1 


2T w /T Ql 


IF, 


Fq2 /T’OI 



(4) 


Equation 4 must be solved by iteration since the mean values of the influence coefficients are functions 
of M 2 , which is not known. The value of M 2 must be assumed, used to compute Fj 0 and Fp and then 
solved for in the finite difference equation. This procedure must be repeated for the length of the duct 
to obtain the distribution of Mach number versus duct length. This is a laborious process for the 
steady-flow case where upstream conditions are usually known. 


In the case for the quasi-steady flows, additional iterations are required which incorporate the 
above procedure as an intermediate step. The only other useful condition besides the initial conditions 
is that, at the exits of the ducts, the Mach number of the flow is unity until the transmitted shock 
reflects from the end of the shock tube and interacts with the drift flow behind the transmitted shock 
wave. (A thorough discussion of the reasoning for this condition is given in Articles 6.3, 6.4, and 8.9 of 
Shapiro (Reference 16).) Therefore, to start an iterative process for determining flow conditions 
before this occurs, a value must be assumed for M 2 , the Mach number of the particle flow between 
the upstream face of the model and the reflected incident shock. (See Figure 2.) The value of M 2 is 
subsonic and provides A 2 /A* where A * is the choked area at the exit plane of the duct and A 2 is the 
upstream channel area. 

The first iteration that must be performed is to provide a value for the Mach number of the re- 
flected wave that would be compatible with the assumed value of M 2 . The expression to compute 
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this value is derived in accordance with Shapiro’s (Reference 1 7) analysis of moving shocks in terms of 
stationary shock formulas: 


_M 2 -{M x ) r 

(■'^3 ) Reflected wave ' U) 

3 C 2 + {M y ) R 

where {M x ) R is the initial Mach number of the reflected wave relative to the oncoming flow. 

C«,>A = U r /C 2 is obtained by iteration. This also provides the stagnation conditions in region 3. 

{A 2 /A*)e = A R /A* ( 6 ) 

where e = plate porosity = A R /A 3 . With the value of A R / A * and the isentropic one-dimensional com- 
pressible flow tables for 7 = 1.4 (Reference 18), all parameters at the duct entry, point R, can be 
obtained. The entry Reynolds number based on duct diameter is then computed. Values of the vis- 
cosity (Reference 19) at the entry conditions are used for this computation. Knowing the Reynolds 
number and the L/d ratio for the duct, an estimate of the average friction factor can be obtained 
(Reference 20) by considering the value at x/d equal to half the incremental length. This estimate is 
valid because the friction factor versus x/d curve (Figure 3) is linear in this initial region. This enables 
the computation of the duct-flow longitudinal temperature profile by the use of Equation 3. When 
the profile is obtained, the successive iterations using Equation 4 are required. If the initial estimate 
of M 3 was valid, M £ should equal unity; if not, the entire procedure must be repeated until con- 
vergence to the desired M R is obtained. 

Limiting Solution Method 

It is readily seen that the aforementioned procedure would be very time consuming and complex 
and that the required additional simplifying assumptions could entail considerable error. An alternate 
method is proposed, which is based on two individual computations of bounding solutions. The 
limiting cases should bracket the physical situation, be considerably less laborious, and prove adequate 
for design purposes. The procedure involves first computing flow through the perforated plate as a 
duct flow for a Rayleigh process of frictionless flow with a T Q change, then computing the duct flow 
for a Fanno process or adiabatic flow with friction. In the case of Rayleigh duct flow, the choked 
section occurs close to the entry at a vena contracta and supersonic flow exists throughout the length 
of the duct. In the case of the Fanno duct, the choked section is at the exit plane and subsonic flow 
occurs for the entire duct length. In both cases, the flow downstream of the perforated plate is super- 
sonic with rapid formation of a base-region flow. When the base-region flows become steady, a 
relatively stable (constant velocity) downstream contact surface with transmitted primary and 
auxiliary traveling waves occurs. The Rayleigh bounding solution is called the “lower limit” because 
its calculations provide lower values of reservoir and transmitted-wave parameters. 

Rayleigh Duct-Flow Formulas (Lower-Limit Solution) 

The following assumptions are made to analyze the flows for this case. The conditions behind 
the incident and reflected waves are based on ideal normal shock theory. Real gas properties are used 
where applicable. The flow behind the reflected wave is compatible with the choked flow at the vena 
contracta near the entrance of the perforated plate. Downstream of the vena contracta, the flow is 
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slightly supersonic and in accordance with the isentropic expansion ratio. The vena contracta is 
assumed to be sufficiently close to the entry plane so that the entire duct can be considered to have 
supersonic conditions. The experiments performed by Weir and others (Reference 21) have shown this 
to be essentially the case. It is assumed that the flow through the duct is frictionless and the duct wall 
is at constant temperature for the brief period of quasi-steady flow. Base region and flows downstream 
of the model are based upon conservation of mass and momentum and compatibility across the down- 
stream contact surface away from the immediate vicinity of the model. 

Upstream Traveling-Wave Formulas 

The following equations for traveling waves are required in order to arrive at stagnation conditions 
for both the Rayleigh duct flow and the Fanno duct flows. They are based on a Galilean transforma- 
tion (Reference 1 7) for converting the traveling plane waves in a stationary coordinate system into 
standing waves in a coordinate system moving at the shock velocity but in the opposite direction. 

Therefore, for the conditions behind the incident wave in region 2 

(M x ), 

M 2 = (M)j (7) 

C 

2 *"'1 

where ( M x ) f = U I /C l is the incident-wave Mach number and (M y )j is the final-state Mach number for 
a normal shock with initial value equal to ( M x ) / . 

_ ) 

Pi P 1 yP X P\ ’ 

Px 



u 2 =M 2 C 2 = M 2 y/yRT 2 g . (9) 

The condition of state ratios in terms of initial and final values, x and y, respectively, are ob- 
tained from the normal shock tables for 7 = 1.4. Equations 7, 8 , and 9 provide all parameters for the 
flow behind the incident shock for the stationary coordinate system or relative to the shock-tube wall. 
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The conditions downstream of the reflected incident wave are obtained from Equation 5, where 
{M x ) R = U R /C 2 . The normal shock relationships hold for the stationary coordinate system as indi- 
cated for region 2 in Equation 8. 


Flow in the Region of the Perforations 


If the Mach number behind the reflected wave is known, the velocity of the reflected wave can be 
obtained from Equation 5; and with the normal shock tables, the steady-state conditions in the reser- 
voir zone can be computed. However, as shown in Figure 4, the drift flow upstream of the model is 
controlled by the vena contracta formed at the square-edged entry. The size of this vena contracta can 
be approximated by considering the empirical work of Koldozie and Van Winkle (Reference 22), who 
expressed an experimental correlation for the discharge coefficients of thick perforated plates in terms 
of hole diameter d, hole pitch P, and a coefficient K. Figure 5 shows the hole pattern specified in these 
terms and also in terms of the parameters considered by Baines and Peterson (Reference 23) and used 
by Dosanjh (Reference 10). Figure 6 is a plot of the correlation constant K versus thickness to hole 
diameter ratio for Reynolds numbers up to 20 000. It can be seen that for L/d values over 2.0 the 
correlation is independent of plate thickness; and, in the Reynolds number range 4000 to 20 000, the 
correlation is also independent of Reynolds number. Since the models of concern have L/d values of 
about 1 0 and Reynolds numbers over 20 000, it can reasonably be assumed that the flat portion of 
the uppermost curve will closely approximate the steady-state performance of the model for the quasi- 
steady-state investigation. (In the absence of data for transient studies of C d , and later heat transfer, 
data from steady-state experiments were used.) The relationship for the family of curves is 


^average 


For the model with a d/P = 0.525 and K = 0.98, C. = 0.919, A„/A* = \/C, = 1.0883, and 



0.10 


A 


3 


A* 


1 A r 

= 4.35327. 

e A* 


( 10 ) 


Both A/ 3 and M R are fixed by the above area ratios and the fluid thermodynamic properties. The 
value of provides the necessary condition to solve Equation 5 for the reflected-wave speed and sets 
the reservoir conditions in region 3. The value of M R provides the entry conditions and the Reynolds 
number based on diameter. Since real fluid properties are being considered here, iterations must be 
performed on M R ,y R ,P R , and T R to compute all conditions at R using one-dimensional compressible- 
flow functions. Tables 7-3, 7-8, and 7-9 of Reference 19 provide additional data required to compute 
Reynolds and Prandtl numbers for the lower Mach number range. For incident-wave Mach numbers 
of 4 and higher, equilibrium thermal conductivity values calculated by Yos (Reference 24) are used 
for obtaining Prandtl numbers. 


To obtain a realistic estimate of the T Q change that occurs in the duct, the following procedure 
is used. Compute the adiabatic wall temperature 
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Figure 4— Details of Rayleigh duct-flow model (typical for each individual 

perforation). 



Figure 5— Model hole pattern after Koldozie (Refer- 
ence 22) and Baines (Reference 23). 



HOLE DIAMETER 


Figure 6— Correlation coefficients for perforated plates 
versus L/d for Reynolds numbers 2000 to 20 000 
(Reference 22). 
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Compute the boundary-layer film temperatures using the static temperature at R for the bulk tempera- 
ture. This should be close because, at L/d =10, the temperature profiles have not been sufficiently 
developed to deviate appreciably from the entrance static values. From Reference 20, shown in Fig- 
ure 7, the value of the parameter X for heat extraction from the gas and 

T,‘X(T w -T r ) + T r , (12) 


where T is the initial temperature of the duct wall, provide a datum reference temperature for evaluat- 
ing the thermal conductivity k ^ of the film at the wall. The value of the mean Nusselt number for L/d 
and the Reynolds number can be obtained from Figure 8. These values enable the determination of an 
effective heat transfer coefficient to be used in computing the convective heat transfer rate into the 
wall by means of Equation 1 3 : 

03 ) 


where A s is the surface area of a single duct. The stagnation temperature change can then be computed 
from the mass flow rate p R u R A R and the specific heat at station R : 
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(14) 


The conditions at the exit plane can now be determined by the Rayleigh line relationships for the value 
of T 0 at E : 

t oe = t or -* t 0 ( 15 ) 


Since the isentropic conditions at R are known from (A R /A *) As=0 , the Rayleigh line unity Mach num- 
ber parametric values can be computed in terms of the isentropic values and the parametric ratios for 
the Rayleigh line at the values of M R . For example, 

(P*W = (Pr Wo ( —] ■ (16) 

\PrJ Ray 

An initial value of M p can be obtained from the stagnation temperature ratio (T 0£ ./7’*) Ray . Iteration 
ofM E ,'y E ,(p E /p*) Ray ,(T E /T*) Ray , and the known value of the stagnation temperature ratio deter- 
mines all conditions at E, the exit plane of the ducts. 


Base-Region Flows 

Using Wick’s relationship (Reference 25) for base pressure to exit-plane pressure ratios in terms 
of area ratio, exit-plane Mach numbers, and Mach number at the point of reattachment of the expanded 
flow, Barthel (Reference 1 1) derived an equation for the base pressure to exit-plane pressure ratio in 
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Figure 7— Variation 


of X with x/d as used in Equation 12 (adapted from Figure 8 of 
Reference 20). 
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Figure 8— Variation of average Nusselt number with x/d and entry Reynolds number (Re ; ) for gas 
flowing in a tube. Uniform wall temperature, uniform initial velocity, and temperature distributions, 
and constant properties are assumed (adapted from Figure 24 of Reference 20). 



terms of the Fanno line impulse functions at the exit plane and the reattached flow plane for the flow 
issuing from the ducts into the full channel. Rearranging the equation and using pertinent steady-state 
experimental base-pressure data presented by Korst (Reference 26) in his work on supersonic axisym- 
metric internal flows (Figure 9), the impulse function at region 4 (the plane of the reattached expanded 
flow downstream of the model) can be obtained: 
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(17) 


where p B fp E = (p B /p QE ) Korst ( p 0E /p E ) Ray • 

Using the value of the impulse function and an estimate of y 4 , the values M 4 ,p 4 , and T 4 can be 
computed and iterated to check against the impulse function value. Initial estimates of p 4 are obtained 
from the equation 
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Initial estimates of T 4 are obtained from the equation 
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Final pressures are computed from the Fanno line pressure ratio, and final temperatures are computed 
using 


'*a\ ! Zl/ffV 
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( 20 ) 


Having the conditions in region 4 and the initial conditions of the shock tube, the transmitted-wave 
patterns can be determined, as shown in the later section, “Downstream Traveling-Wave Formulas.” 


Fanno Duct-Flow Formulas (Upper-Limit Solution) 

Since the computation of the duct flow involving real gas properties requires extensive iteration, 
the Fanno line solution is based on the gas having a constant 7 = 1.4, and table values for one-dimen- 
sional isentropic, normal shock, and Fanno line functions are used (References 18 and 27). For this 
case, the flow is assumed to be adiabatic and friction is considered in order to determine the entry 
Mach number M R that will produce choking at the exit plane E. To start the computation, an initial 
value of the duct entry Reynolds number must be assumed in order to establish an average value of the 
friction factor. This provides a value of 4 fL/d from which the Mach number at R can be obtained 
using table values of the one-dimensional compressible flow functions for a Fanno line. For the same 
value of the entry Mach number and the isentropic area ratio corresponding to this number, the isen- 
tropic Mach number for the particle flow in region 3 (the stagnation region) can be obtained since 
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Figure 9— Experimental base-pressure ratios for super- 
sonic axisymmetric internal flow (adapted from Figure 27 
of Reference 26). 
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As=0 


The value of M 3 corresponding to the resulting 
area ratio is then fixed for use in the iteration to 
obtain the reflected-wave speed using Equation 5. 
As in the case of the Rayleigh duct flow, this 
establishes the reservoir conditions upstream of 
the entry and, using data from tables for isen- 
tropic, subsonic conditions, all parameters at R 
can be obtained. The Reynolds number at R can 
then be computed and checked against the 
initially assumed value; if it is such that the 
average friction factors do not change signifi- 
cantly, the original assumption is acceptable. 
Otherwise, the iteration process must be repeated 
for the new value of the Reynolds number. 

Using Fanno line functions and M E = 1 at 
the exit plane of the model, all parameters at E 
can be determined. To determine the conditions 
in region 4, Kochendorfer’s analysis (Reference 
14) for a sonic jet expanding into an annular 
enclosure that is closed at the upstream end and 
open at the downstream end was used. This 
analysis assumes a steady adiabatic flow between 
the sonic exit of the duct and the station at which 
the sonic jet expands and attaches to the annular 
enclosure and is directly analogous to the base 
flows expected during the relatively steady 
period of the Fanno expansion of the Fanno 


line solution. From the equations for conservation of mass and momentum, an expression for the 
Mach number at the plane of reattachment in terms of the porosity ratio of the plate is obtained: 



MU - 7 Ml , ( 22 ) 


and the pressure in the same plane in terms of the exit-plane stagnation pressure is obtained: 


p 4 ey/[ 2 /(y+iy]<rHMy-i) 


Poe M 4 y/l + [( 7 - l)/2]M 2 4 


( 23 ) 
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This analysis provides the lower limit of the pressure ratio required for the flow to fill the 
channel and assumes for steady-flow conditions a normal shock at the plane of reattachment. It does 
not consider entropy losses due to expansion or boundary-layer effects. However, since a thick turbu- 
lent boundary layer is assumed to exist and in essence forms a sonic nozzle at the exit plane, the area 
ratio can be based on exit plane to free channel area or AjA* instead of the geometric area ratio; this 
provides a slightly higher value of Mach number in region 4. The compatibility of region 4 with the 
downstream wave patterns can then be determined. For convenience, a plot of area ratio versus 
region 4 Mach number based on the Kochendorfer analysis is shown in Figure 1 0. 

Downstream Traveling-Wave Formulas 

As indicated earlier, there are three possible auxiliary-wave patterns that can be associated with 
the transmitted shock. The formulas indicated in the following sections will provide the means of 
establishing compatibility at the transmitted contact surface of conditions in the associated regions. 

Auxiliary Left-Facing Standing Shock Wave 

At lower incident-wave Mach numbers (on the order of 2 or less), the possibility of an auxil- 
iary left-facing wave should be considered. For this case, the wave should form at or slightly 
downstream of the point where the base-region boundary stream perimeters intersect and refract. 
The short region of refraction forms a recompression zone and the complex three-dimensional field 
rapidly coalesces to form a one-dimensional standing wave. The particle flows on either side travel 
downstream, and the elevated pressures and temperatures are determined by the normal shock 
relationships. The upstream Mach number is ( M x ) A = M 4 , and the Mach number downstream of 
the standing wave is ( M y ) A - . Therefore all conditions at region 5 can be obtained from 



MACH NUMBER M. 

4 

Figure 10— Expansion ratio versus Mach number according to Kochen- 
dorfer's analysis (Reference 14). 
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region 4. Across the contact surface, for quasi-steady conditions, p $ = p 6 and u $ = u 6 . These are the 
compatibility constraints that must be satisfied in determining the speed of the transmitted waves. An 
iterative process is required to arrive at p 6 and u 6 working back from the region of initial conditions, 
region 1. If a compatible solution does not exist for the standing wave, a balance will be achieved for 
some value of that is different from that determined from base-region considerations. 

The region 6 iteration is performed by assuming the Mach number of the transmitted wave 
(M x ) t = U T /C X . The normal shock relationships provide (M y ) T> p 6 = p y9 and T 6 = T y since p x = p x 
and T x = T x . By transforming the coordinate system from a stationary one to one moving at the 
transmitted shock velocity, an expression for particle velocity u 6 is obtained: 

u 6 =(M x ) T C l -(M y ) T C 6 . (24) 

The existence of the auxiliary standing wave can then be tested. 


Auxiliary Left-Facing Traveling Shock Wave 


The left-facing auxiliary traveling wave is the most likely to satisfy the transmitted contact surface 
compatibility constraints for this incident-wave Mach number range, and its existence has been experi- 
mentally verified by several investigators. The wave, although left-facing, will travel downstream with 
reference to a stationary coordinate system at the velocity U A - u 4 and its rate will increase with an 
increase in incident-wave strength. 

The value of M $ in this case is obtained by assuming the value of (M X ) A = U A /C 4 and using the 
relationship 


M a - (M ) 


m 5 = 


x'A 


5 C 4 


+ (M y ) A ■ 


(25) 


Using the normal shock relationships, p $ = y p x p A , T s = y T x T A , 5 C 4 = y/yT x , and u $ =M $ C S . Since 
p 5 = p 6> pJP\ determines the Mach number of the transmitted wave. By means of the normal shock 
relationships and the initial conditions, all parameters in region 6 can be computed. The particle 
velocity in this region, u 6 , can then be computed using Equation 24. If u 6 is larger than u $ , a smaller 
value of (M x ) a is assumed and vice versa until u $ checks u 6 within 1 percent. 


Auxiliary Left-Facing Rarefaction Wave 

Barthel (Reference 1 1) suggested the possibility of the formation of a transmitted auxiliary rare- 
faction wave as incident-wave strengths are increased. His reasoning was based on the gradual lowering 
of velocity of the left-facing auxiliary shock as incident-wave strength increased. Assuming the velocity 
was sufficiently lowered so that the auxiliary shock strength approached a sound wave, this possibility 
could follow. However, he did not experience this in the low Mach number range that was tested. It 
is doubtful that this will occur even at higher incident-wave Mach numbers, since analysis predicts 
appreciably higher wave strengths for a rarefaction wave then it does for an auxiliary traveling shock 
wave which should not occur with sound waves. 
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This analysis was based on relationships between sound speed and pressure at a rarefaction wave 
and on maintaining compatibility at the transmitted contact surface. Again an iterative process is 
required to utilize the fundamental equations. A value of p s must be assumed to compute the sound 
speed in region 5. Since motion of sound waves is isentropic, the Mach number in the region can be 
computed using the stagnation pressure of region 4. The sound speed and Mach number provide the 
particle velocity in the region, which must be balanced against the particle velocity of region 6 as com- 
puted in the previous section. The equations required are as follows: 



(26) 
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( 28 ) 


The pressure in region 5 must balance the pressure in region 6. Therefore, the primary transmitted- 
wave shock strength is specified and all parameters in region 6 can be computed based on the initial 
conditions of region 1. The particle velocity n 6 computed from Equation 24 must equal u 5 or the 
iteration must be carried further. 


The numerical results from the computations described in the preceding sections are presented 
and discussed later. 


APPARATUS AND EXPERIMENTAL PROCEDURE 

To assess the correlation of the theoretical calculations with actual flow behavior, an experimen- 
tal program was conducted. A shock tube was designed to generate plane traveling waves in the range 
of incident-wave strengths previously discussed. Particular attention was paid to surface finish, sur- 
face discontinuities, rigidity, and instrumentation capability. This section describes in detail the basic 
facility and its instrumentation, the models used, and the test procedure. 

Shock-Tube Design 

The shock tube. Figure 1 1 , is a high-strength, thick-walled, rectangular steel box of uniform 
cross section having inside dimensions of 1.50 in. high by 2.00 in. wide. The first-stage helium driver 
section consists of two 24-in. -long flanged units bolted to make up a 48-in. -long reservoir. If required, 
a third 24-in. -long unit can be inserted to provide a 72-in. driver. The second-stage driver is a flange- 
less unit 2 in. long and is clamped by four massive C-bars between the first-stage driver and the driven 
section. Two identical pure aluminum diaphragms are interposed between the 2-in. section and the 
clamping flanges. The nitrogen-filled driven section consists of five 24-in. -long flanged units. The 
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Figure 1 1— Overall view of shock-tube laboratory. 


model injector plates are clamped between the third and fourth units, which have lO-in.-long full- 
height glass or plastic windows adjacent to the model faces to enable optical recording of interaction 
phenomena. Since the incident wave can travel 72 in. prior to impinging on the face of the model, 
it should have ample time to develop a uniform strength across the tube cross section. The proportion- 
ing of the basic tube and the detail design was in accordance with generally accepted practice 
(References 28 to 30). 

The double diaphragm arrangement is used in conjunction with cruciform-milled aluminum plates 
to eliminate the need for diaphragm piercers. The 2-in. section is pumped to about one-half the pres- 
sure of the first-stage driver. When the pressure of the 2-in. section is exhausted to atmospheric 
pressure by means of a solenoid valve, the upstream diaphragm ruptures and dumps the full pressure 
on the downstream diaphragm. The downstream diaphragm then rapidly ruptures with minimal 
deformation because of the higher loading rate. The diaphragms petal cleanly without tearing and the 
repeatability of the breaking load can be kept within 5 percent on properly fabricated units. 

All major shock-tube parts are made of heat-treated SAE 41 50 steel plates that were machined 
and ground to size and have at least an 8-pin. rms internal finish. The plates are electrodeless, nickel 
plated, and electroblackened to provide a hard corrosion-resistant dull coating on the surfaces. The 
plates and flanges were assembled with high-strength bolts and dowels. The assembled sections are 
supported on individual carriages which permit lateral positioning. The carriages roll on hardened and 
ground guide rods by means of chain-linked roller bearings. The guide rods rest in V-grooves milled 
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into the rigid structural steel base; the base is mounted on heavy-duty swivel casters and jacks for 
portability and positioning. Instrument ports are located in the top and bottom plates of all sections. 
They are Vi in. in diameter and are filled by O-ring sealed plugs mounted flush with the inside wall. All 
temperature and pressure measuring instrumentation is inserted into the system by means of these ports. 

The upstream and downstream test sections are similar in construction to the other units except 
that thick fused-silica or cast acrylic windows (optically flat) provide a clear view of the interior of the 
flow channel 10 in. in either direction from the faces of the model when it is clamped between the two 
test sections. There are instrument ports every 2 in. starting at the model face for the first foot of 
length in the top and bottom plates. These provide a flexible arrangement for tranducer location and 
for triggering auxiliary electronic circuitry. Fiducial wire grids clamped outside the glass windows 
provide a convenient scale of reference for optical determinations of wave speeds. Massive T-bars, 
shown in Figure 12, restrain the cantilevered ends of the upper and lower plates at the model faces. 

This is necessary in order to provide the support lost by having a split clamp flange at the model faces; 
the flange was split to offer an unobstructed view of the model-gas interface. 

The shock tube is manifolded to vacuum pumps and to nitrogen and helium supplies by means of 
a graphic control panel. The panel controls enable the evacuation of all three sections of the shock 
tube, pressurization of the driver sections using either an 800- or a 3500-psi helium supply, and 
pressurization of the driven section using a low-pressure (80-psi maximum) nitrogen supply. The panel 



Figure 12— View of test section showing split flanges, T-bars, model, fiducial grids, and instrumentation. 
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contains master gages for indicating vacuum and pressures for system scavenging and charging purposes. 
The gages and supply lines are piped to permit in-line calibration with primary and secondary standards. 

The gases are supplied from bottle banks and cleaned by passing through purifier units prior to 
regulation. The higher pressure helium is obtained by boosting gas regulated at 800 psi to 3500 psi in 
a diaphragm compressor and storing in high-pressure receivers (1772 in. 3 ). The shock tube can be 
pumped down to 0. 1-mm-Hg pressure in less than 10 min using a small two-stage vented-exhaust 
vacuum pump rated at 140 liters/min. For higher pumping rates, it can be manifolded to two tandem- 
arranged vacuum pumps each having a nominal rating of 1600 ft 3 /min. 

Instrumentation 

The instrumentation used with the facility can be categorized as either steady-state or transient 
measuring systems. The steady-state instruments include those required to establish initial conditions 
in the shock tube. The transient measuring systems are those that provide the data required to assess 
interaction phenomena produced when the generated traveling waves impinge on the model faces. 

Steady-State Measurements 

The steady-state measurements of pressure in the various sections of the shock tube are made with 
solid-front, temperature-compensated units. The vacuum that scavenges the system prior to each run 
is measured by an S^-in. 0- to 760-mm-Hg (1-mm-Hg divisions) gage using ambient pressure as the 
datum point. This same gage is also used to measure the low pressure to which the nitrogen-filled 
driven section is charged. To monitor the driven section pressure when the more precise vacuum gage 
is valved off, a 6-in. 0- to 1000-psi ( 1 -psi divisions) absolute pressure gage is connected to the driven 
section supply line. 

Two 8!/2-in. gages are manifolded to each driver section to cover the complete range of charging 
pressures. The short second-stage driver is measured by 0- to 250-psi (0.5-psi divisions) and 0- to 
2500-psi (5-psi divisions) units; the long first-stage driver is measured by 0- to 500-psi ( 1 -psi divisions) 
and 0- to 5000-psi (10-psi divisions) units. High-pressure hand valves enable the selection of the gages 
suitable for the test pressure levels. 

The initial temperature of the gas at rest is assumed to be the same as that temperature monitored 
within the channel wall at a distance of 0.061 in. from the interior surface. This assumption is reason- 
able since adequate care is taken to insure that the equilibrium wall temperatures are close to ambient 
conditions; the small volumes of gas within the channel should quickly stabilize to the wall tempera- 
tures. Type K (chromel-alumel) welded thermocouples are located in four instrument plugs at 
stations 37, 60, 1 16, and 10D of the tube. Two additional thermocouples are inserted into the model. 
The thermocouple extension leads are terminated at a thermoelectric manual-switching unit which 
enables the selection of a single output for display on a digital thermocouple indicator. When all 
temperature channels agree within 1 F°, the channel is considered to be at thermal equilibrium and 
ready for a test run. Figure 13 is a view of the control panel and the digital temperature-indicating 
system. 
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Figure 13-Shock-tube pressure control panel and digital temperature- 

indicating system. 


Transient-State Measurements 

Three major modes of data acquisition are utilized to measure the flow fields prior to and during 
interaction with the models. They are 

( 1 ) High-responsc pressure transducers to acquire pressure histories 

(2) High-speed electronic counters to measure average wave speeds at discreet positions in 
the tube 

(3) High-speed framing camera to provide a wave-speed history for a distance of 10 in. upstream 
and downstream of the model and to photographically record interaction phenomena 
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High-Response Pressure-Measuring System 

The components of the pressure-measuring system include quartz pressure transducers having a 
nominal sensitivity of 0.35 pCb/psi and a resonant frequency greater than 400 000 Hz. The transducer 
is flush mounted into the shock-tube wall by means of a brass instrument port plug machined to accept 
the transducer, its metallic seal, and its connector adaptor. The brass plug is retained by a 5/8-18 UNF 
thread at the outer wall surface, but is sealed at the inner wall by an O-ring on its 0.500-in. diameter. 
The charge signals generated by the shock pressure steps are converted into low-impedance voltage 
signals compatible with cathode-ray oscilloscopes by means of electrostatic charge amplifiers. Two 
channels upstream and two channels downstream of the model are used to record pressure histories. 
The oscilloscopes are wide-range dual-beam units providing a frequency capability in the DC to 30- 
MHz range. They permit individual or synchronized control of the time base plug-in units; the vertical 
channels utilize separate plug-in preamplifiers. The upper beam preamplifier is a fast-rise, high-gain 
unit with a deflection capability of 5 mV/ cm to 20 V/cm having a passband of DC to 30 MHz and a 
rise time of 12 ns; a square wave is fed into this channel by a pulse generator to provide a voltage and 
frequency calibration signal for data reduction of the lower-beam pressure traces. The lower-beam 
plug-in preamplifier is a four-trace unit having a DC to 20-MHz passband and an 1 8-ns rise time, with 
a deflection factor of 0.02 to 10 V/cm. The two pressure transducers on the upstream side of the 
model use the first two channels of this preamplifier in the chopped mode of operation. An identical 
arrangement exists for the two pressure channels on the downstream side of the model with a second 
oscilloscope whose upper beam is also linked to the calibration pulse generator. 

A total of nine pressure transducers is distributed along the length of the shock tube. In addition 
to the four transducers that provide pressure histories and trigger complementary instrumentation, 
five other units are used for triggering purposes only. Figure 14 shows the transducer installations and 
the close-coupled charge amplifier racks. Figure 15 schematically indicates the overall instrumentation 
arrangement. 



Figure 14— Pressure transducer ana charge amplifier installation for data acquisition and 

triggering functions. 
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Wave-Speed Measurement 


The average wave speeds are measured in selected locations upstream and downstream from the 
model. The incident-wave speed is averaged over a 10-in. interval 1 2 in. and 2 in. upstream of the 
model. The reflected-wave speed is measured over a 2-in. interval 2 in. upstream of the model. The 
transmitted-wave speeds are measured over an 8-in. interval 2 in. downstream from the model and 
over a 1 0-in. interval 1 2 in. downstream from the model. 

The time intervals for traversing these distances are recorded on a multichannel counter. This 
counter system has two master clocks, a common master start, individual stops, individual resets, and 
master reset capability. The clock frequency is 10 MHz with a minimum timing capability of 0.3 /is. 
The counters are started when the incident wave excites the transducer at station 90, which is 32 in. 
upstream of the model. The transducer signal is amplified by the charge amplifier and then fed 
directly to the start input of the counter bank. The charge amplifier range is set to produce about a 
1-V signal at each incident-wave speed for reliable triggering. The stop circuits have additional 
auxiliary electronics since the other charge amplifiers are set at only one range and are not changed. 
This is done to avoid a variable response due to the effect of amplifier feedback capacitance on unit 
time constant and to match amplifier delays. To provide the additional amplification for triggering 
on these channels, the outputs of the charge amplifiers are fed into low-gain, high-response operational 
amplifiers Qu « 20) and the outputs of these are then used to trigger the counter channels that 
monitor downstream traveling-wave speeds. To measure the reflected-wave speeds (upstream travel- 
ing), the outputs of the operational amplifiers at stations 1 18 and 120 (4 in. and 2 in., respectively, 
upstream of the model) are fed to pulse-counter circuits that provide a fixed output to the counter 
gates when the second-higher-level pressure pulse of the reflected wave triggers these units. Figure 16 
shows the rack containing the counter system, the triggering and gating unit (the pulse counters), 
and the conditioner unit (operational amplifiers). 

Having direct digital indication of incremental elapsed time from the counter bank and the 
incremental distances between transducer stations, the average wave speeds can be readily determined. 

Optical Instrumentation 

The optical recording system consists of an//8 folded Z-type schlieren system, two pairs of 
optical flats in the vertical walls of the shock-tube test section, appropriate light sources, and a 
continuous writing streak and framing camera. The schlieren components are supported on two 
massive (1400-lb) pedestals that are aligned parallel to the shock tube. They include off-axis parab- 
olas and produce a fully collimated lO-in.-diameter working beam that is perpendicular to the shock- 
tube longitudinal axis. The shock tube is illuminated by arc lamps in the vicinity of the optically 
flat windows of the test sections. Since the flow channel is mounted on movable carriages, it can be 
positioned so that the upstream window, injector plate model, or downstream window are centered 
in the working beam. The 1.5-in.-high by lO-in.-long glass windows and mirror system components 
are diffraction limited (ground and polished to within one-tenth the wavelength of green light, 

5400A). The off-axis angle of the fused silica parabolas is less than 7°. The steady-state light source 
used for alignment purposes is a mercury-vapor unit. High-pressure 80-W short-arc xenon flash lamps 
are used for short-duration illumination exposures for the camera. They are energized by inductance- 
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Figure 16— Wave-speed measuring system consisting of chronograph units, 
operational amplifiers, and pulse-counting gating unit. 


capacitance pulse-forming networks, similar to those described by Mauldin and Compton (Reference 
31), to produce square wave pulses of light having durations of 100, 400, and 1200 ps. 

The camera system, modified to run at 1 2 000 to 2 000 000 frames per second, is used in con- 
junction with an auxiliary field lens assembly (Figure 1 7) so that the field of view is about 10 in. in 
diameter to match the working beam size. The 35-mm format provides a nominal frame size of 16 by 
25 mm. With the 360-mm-focaI-length objective lens set at infinity, the 10-in. field should com- 
pletely fill the frame. The rotating mirror turbine is driven at a rate compatible with the selected 
light pulse duration to prevent overwriting on the two 40-frame film strips that comprise the record 
of a single run. 

The light source is triggered by the output of the pressure transducer located 12 in. upstream 
from the model (station 1 10). The transducer charge is converted to a voltage signal in the charge 
amplifier and then fed to the high-voltage trigger unit which provides a 200-V maximum negative 
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Figure 17— The continuous streak and framing camera system (camera turbine and shutter controls shown 

at inset). 

step output. The negative step triggers the first channel of a two-channel delay generator unit. The 
delay generator output is a positive pulse of 250 V; it is used to trigger the light-source power supply 
and the camera turbine air cutoff control. The output of the high-voltage trigger unit is also fed into 
the second channel of the delay generator console. The output of this second channel triggers the 
square wave pulse generator (frequency and voltage calibration signal) that is fed to the oscilloscope 
preamplifiers, which in turn triggers the sweeps of the two oscilloscopes that record the pressure 
traces on positive photographs. 

The 35-mm film strips provide data on wave speeds, wave acceleration, and a photographic 
record of the flow-field development. When examined with the pressure traces, they provide a com- ’ 
plete history of the transient interaction phenomena. 

Kodak Plus-X Pan 35-mm film was used for all pictures and the developer was Diafine dual bath 
with a developing cycle of 10 min in bath A and 6 min in bath B at 70° F. 

Test Configuration 

As indicated in the discussion of the shock-tube design, the models are in the form of plates that 
are clamped between the split flanges of the upstream and downstream test sections. An 0.030-in.- 
thick rubber gasket is seated around the 0.0 80-in. -high by 1.50-in. by 2.00-in. boss on each model 
face. This provides a vacuum-tight seal between the models and the test section faces when properly 
assembled and bolted. Dowels inserted into the ends of the horizontal plates of the test sections 
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provide alignment of the models and the shock tube. The models, shown in Figure 18, are machined 
from aluminum alloy 2024-T4. The hole patterns were produced by precision drilling, reaming, and 
lapping to size. Edges were maintained square and sharp by grinding and lapping the faces of the 
plates. 

Two models were considered for the experimental program. Both have nominal Lid values of 10 
and porosity factors of 0.25. The hole sizes were selected to provide a hole area scaling factor of 10 
to 1 between models. One model has 179 effective passages of 0.078 in. in diameter by 0.875 in. long 
resulting in an L/d value of 11.15; the other model has 1 7 effective holes 0.250 in. in diameter by 
2.500 in. long with L/d equal to exactly 10. Thermocouple wells drilled from the top surface enable 
measurement of the model temperatures on the vertical centerline to within 0.060 in. of the upper- 
most flow passage. The wells are located 0.187 in. from the model faces on the 0.875-in.-thick 
model and 0.500 in. from the faces of the 2.500-in.-thick model. Number 30 AWG chromel-alumel 
thermocouples mounted on slotted 0.125-in.-diameter Micarta rods are fitted to the wells so that their 
hot junctions make intimate contact with the flat bottoms of the wells. 

Test Programs 

For each model, experiments were run for selected positions of the model relative to the schlieren 
beam with the incident-wave Mach number varied over the nominal range of 2 to 4.6. To test over 
this wide range of Mach number with maximum safety, cast acrylic windows were used instead of the 
glass. Since it was not feasible to polish the plastic to the flatness of glass, there was a resulting loss of 
resolution at the faces of the model due to edge diffraction. 

The initial position presented the flow field of the incident and reflected waves at the upstream 
model face; the second position centered the model in the schlieren beam to show upstream and 
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downstream fields in close proximity to the model; the third position permitted examination of the 
transmitted-wave patterns downstream of the model face. A minimum of four values of incident- 
wave Mach number were tested in the designated nominal range. At each value of Mach number, five 
repetitive tests were performed to assess variability. Schlieren beam position was varied to obtain 
optimal data. In general, initial conditions, counter-wave-speed data, pressure transducer output, and 
high-speed photographic records of the various flow fields were acquired for each of the test runs in 
the ranges Mach 2 to 4.6 for the 0.078-in.-hole-diameter and Mach 2 to 3.75 for the 0.250-in.-hole- 
diameter model. Thus for a full range of tests for each model, 20 to 25 runs were performed. 

Experimental Procedure 

The following is a description of the sequential operations that were used to perform a single test 
run with the 1.5- by 2.0-in. shock tube. 

The selected model was clamped between the test sections and a leak test was performed to 
guarantee that the assembly was within tolerance; i.e., 5-mm-Hg maximum pressure rise in 15 min 
after evacuation to less than 1-mm-Hg absolute pressure. 

The amplitude and frequency of the pulse generator output was adjusted to selected values using 
secondary standards. The oscilloscope time base and amplification ranges were set using the calibrated 
output of the pulse generator. Counter attenuations were set and counter “start” trigger levels were 
adjusted. 

Shock-tube and model temperatures were permitted to stabilize to within 1 F° of ambient and 
recorded. Vacuum and pressure gages were zeroed. The appropriately sized diaphragms were seated 
at both faces of the 2-in. driver section with stress grooves facing downstream, and the assembly was 
clamped into the shock tube using the C-bars (Figure 19). Bolt torques were kept lower than 90 per- 
cent of the maximum allowed. The driver and driven sections were evacuated until the vacuum gage 
varied from the ambient barometric pressure by less than 0.5 mm Hg. The driven section was then 
charged with dry, filtered nitrogen (dewpoint, -98° F) to a pressure of 5 psi absolute. The two driver 
sections were brought up to the pressure required for the 2-in.-long section using dry, filtered helium. 
Finally, the first-stage driver was brought to full pressure. 

The framing camera was loaded with film strip from the cassette loader magazines; the camera 
turbine speed was adjusted to correspond to the duration of the selected arc lamp flash. The trans- 
ducer charge amplifiers were grounded, oscilloscope sweeps reset, counters reset, and delay generators 
reset. The room lights were turned off and shutters of the framing and oscilloscope cameras opened. 
The solenoid dump-valve push-button switch was operated. The scope camera shutters were closed, 
the room lights were then turned on, and the framing camera strips retrieved into the cassettes for 
processing. Counter readings were recorded. Oscilloscope records were removed and properly 
annotated with run number, incident-wave Mach number, model number, beam position, and pulse 
generator input wave amplitude and frequency. 

DISCUSSION OF RESULTS 

The following discussion deals with the results of the theoretical consideration of the Rayleigh 
and Fanno flow models in terms of the compatibility of the three possible auxiliary-wave patterns 
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F igu re 1 9— Secondary driver section showi ng vent valve and C-bar clamps. 


that can be associated with the transmitted shock. The three solutions corresponding to these wave 
patterns are referred to as indicated earlier: ( 1 ) the auxiliary left-facing standing shock, (2) the 
auxiliary left-facing traveling shock, and (3) the auxiliary left-facing rarefaction wave. The experi- 
mental data required to validate the theoretical results, qualitatively and quantitatively, are presented 
and a comparison is made. Finally, an error analysis of the data systems used in the experimental 
program is detailed to enable a realistic assessment of the experimental verification of the theory. 

Pressure ratios and wave speeds are the parameters used in the discussion of the results. These 
were selected primarily because they were major parameters computed in the theoretical portion of 
this work and correlation could be directly obtained from the experimental data. Parameters such 
as Mach number were computed from these primary measurements. Although temperatures were 
computed for the theoretical analysis, they are not presented in the discussion because temperature 
histories were not measured in the experimental program. 
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Theoretical Results 


Two sets of theoretical computations for the primary model configuration (0.07 8-in. -diameter 
holes) were performed for the range of incident-wave Mach numbers previously discussed; one set was 
in accordance with the Rayleigh (or lower limit) duct-flow model and the second set was according to 
the Fanno (or upper limit) duct-flow model. In addition, a set of computations for incident-wave 
Mach number 3 was performed for the secondary model configuration (0.250-in. -diameter holes) in 
accordance with the Rayleigh model to examine the influence of a higher Reynolds number on the 
results. In all computations, the compatibility of the three possible types of auxiliary transmitted 
waves was examined. 

The results of the analysis with the Rayleigh duct-flow model are presented in terms of computed 
quasi-steady pressure ratios in the flow regions upstream of the test model, at the model faces, and 
downstream of the model in Table 1 . 

For the range of incident-wave Mach numbers considered, 2 to 5, these computations did not 
indicate the validity of the left-facing auxiliary standing wave downstream of region 4. The solutions 
for the auxiliary standing-wave cases shown in the table may represent an interim state during the 
transient period of base-region pumping that could exist at a Mach number lower than that which must 
theoretically exist in region 4 to satisfy continuity and momentum. For these cases, at the conclusion 
of the initial base-region transient, the quasi-steady stable auxiliary solution would be more aptly 
represented by one of the second type (left-facing traveling wave) or one of the third type (left-facing 
rarefaction wave). 


Table 1— Theoretical pressure ratios versus incident-wave Mach number for Rayleigh duct-flow model. 

Mach 

Number 2Pl 3P i rPi 

Primary 2 4.50 12.54 4.89 

(p j = 5 psia and 
T l = 535° R) 

3 10.33 40.19 13.72 

4 18.50 87.41 31.59 


5 29.00 150.60 54.70 

Secondary 3 10.33 42.12 15.02 

(Pj = 5 psia and 
T x = 535° R) 


eP i 

bP i 

4^1 

Solution 

6 Pi 

3.67 

0.39 

0.63 

i 

2.82 




2 

3.05 




3 

4.04 

10.27 

1.13 

1.77 

1 

4.25 




2 

6.27 




3 

8.20 

21.61 

2.20 

4.08 

1 

6.14 




2 

11.43 




3 

15.60 

42.70 

4.56 

7.62 

1 

9.62 




2 

17.62 




3 

21.40 

11.98 

1.30 

2.02 

1 

4.48 




2 

6.74 




3 

8.72 
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To arrive at these interim solutions of the first type, the conditions of state in region 4 were 
maintained and the Mach number was assumed to take on a lower value that would satisfy compati- 
bility across the transmitted-wave contact surface. It should be noted that for an incident-wave Mach 
number of 2, the region 4 Mach number had to be lowered less than 7 percent from the value obtained 
from the impulse function (Equation 1 7) to enable this hypothetical compatibility; at an incident- 
wave Mach number of 5, the region 4 Mach number had to be lowered by 46 percent to produce 
similar agreement. In essence, at the low end of the Mach number range, the Rayleigh theory predicts 
the possibility of quasi-steady compatibility of the first type; but, as the strength of the incident wave 
increases, the theory definitely indicates that either the traveling wave or the rarefaction waves are 
the theoretically compatible auxiliary phenomena. 

It should also be noted that the solution of the third type results in highest final pressure ratios, 
so that as incident-wave strengths are increased and the strength of the auxiliary wave decreases, the 
compatibility criteria, based on one-dimensionality of the quasi-steady fields, indicate that the 
appearance of a rarefaction wave (as predicted by Barthel (Reference 1 1)) is not likely since the 
pressure ratios predicted for this solution are significantly higher than that required for the left-facing 
auxiliary traveling shock wave. Since rarefaction waves are sound waves with relatively small pressure 
perturbations across their fronts, it is doubtful that this phenomenon would be evidenced as a 
simple auxiliary wave except at much higher incident-wave strengths than are being considered here. 

The set of computations for the secondary test model demonstrated that the greater Reynolds 
number was no advantage as far as heat transfer was concerned. The higher Nusselt number was 
negated by larger hole diameter so that the net effective heat transfer coefficient was 14.5 percent 
lower, thus resulting in a lesser drop in stagnation temperature. This is reflected in the pressure ratios 
shown in Table 1 . The duct exit-plane pressures are 16.66 percent higher for the model containing 
the 0.250-in. holes; transmitted pressure ratios for the three types of solutions are 5.4, 6.4, and 7.5 
percent higher. 

Table 2 indicates the velocities of the various waves as a function of the incident-wave Mach 
number. The final column of the table indicates that although the second solution wave speeds in- 
crease negatively, the particle velocities in region 4 increase at an even greater rate; thus the auxiliary 
waves (when referenced to a stationary coordinate system) appear to move downstream at a monoton- 
ically increasing rate. 

The theoretical pressure ratios for the Fanno duct-flow model are shown in Table 3. They are, 
in general, significantly higher than those for the Rayleigh computations. 

Theoretical wave speeds, shown in Table 4, in general, are slightly higher except in the case of 
second solution auxiliary-wave speeds. These are of comparable value for all incident-wave Mach 
numbers except 2. At this value, the second solution wave speed becomes zero or reverts to a stand- 
ing wave. This is reinforced by the values of the primary transmitted-wave pressure ratios of the first 
and second solutions; their values are less than 1.4 percent apart. In addition, the solution of the 
first type indicates the full compatibility for a left-facing standing wave as an auxiliary solution in 
terms of all conditions of state and velocity at region 4. 

Therefore, in accordance with the Fanno method of computation, at the initial incident-wave 
Mach number, the auxiliary transmitted wave (after the brief initial transient period) should be the 
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Table 2— Theoretical wave speeds versus incident-wave Mach number for Rayleigh duct-flow model 

(stationary coordinate system). 



Mach 

Number 

Incident 

Reflected 

Auxiliary- 

Transmitted 

Auxiliary 

Model 

U > 

U R ~ u 2 

Wave 

u T 

U A ~ W 4 


(fps) 

(fps) 

Solution 

(fps) 

(fps) 

Primary 

2 

2306 

-943 

i 

1845 


{p { = 5 psia and 




2 

1916 

90.7 

7\ = 535° R) 




3 

2190 


3 

3459 

-1040 

1 

2242 






2 

2708 

579 





3 

3099 



4 

4612 

-1171 

1 

2681 






2 

3663 

1325 





3 

4266 



5 

5765 

-1368 

1 

3363 






2 

4535 

1542 





3 

5032 


Secondary 

3 

3459 

-1042 

1 

2301 


(p x = 5 psia and 




2 

2802 

619 

T { = 535° R) 




3 

3102 



Table 3— Theoretical pressure ratios versus incident-wave Mach number for Fanno duct-flow model 

for the primary test model. 


Mach 

Number 

lPi 

3 P 1 

rP l 

eP 1 

aP\ 

Solution 

f>P \ 

2 

4.50 

12.66 

9.72 

5.90 

0.59 

1 

2.94 







2 

2.90 







3 

3.74 

3 

10.34 

43.92 

33.18 

20.84 

2.10 

1 

4.66 







2 

7.16 







3 

9.20 

4 

18.50 

94.38 

71.10 

44.90 

4.50 

1 

7.20 







2 

12.86 







3 

16.34 

5 

29.00 

162.78 

122.64 

77.44 

7.76 

1 

10.88 







2 

20.74 







3 

26.10 
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Table 4— Theoretical wave speeds versus incident-wave Mach number for Fanno duct-flow model for 

the primary test model (stationary coordinate system). 


Mach 

Number 

Incident 

(fps) 

Reflected 
U R - u 2 
(fps) 

2 

2306 

-954 

3 

3459 

-1108 

4 

4612 

-1325 

5 

5765 

-1569 


Auxiliary- 

Transmitted 

Auxiliary 

Wave 

u T 

U A ~ M 4 

Solution 

(fps) 

(fps) 

1 

1882 


2 

1872 

0 

3 

2102 


1 

2345 


2 

2888 

613 

3 

3267 


1 

2897 


2 

3853 

1105 

3 

4337 


1 

3548 


2 

4881 

1531 

3 

5471 



left-facing standing wave; for higher values of the incident-wave strength, the auxiliary transmitted 
wave should be the left-facing traveling wave. 

Hence, both the Rayleigh and Fanno theoretical flow models tend to predict similar auxiliary- 
wave patterns and parametric trends for the same incident-wave Mach number. 

It should be noted that, in the Fanno limit series of computations, the area ratio used in Kochendor- 
fer’s method to compute the Mach number in region 4 was not the geometric ratio based on AJA E , 
but instead, since the flow was assumed to choke at the exit plane, A 4 /A * was used. This is closer to 
a real case than the geometric area ratio since one would expect a significant thickness of the boundary 
layer at the exit plane. The progressive growth of the boundary layer along the length of the duct 
contributes to the choking of the duct and is the controlling factor in adjusting the mass flow and 
upstream stagnation conditions in region 3. A comparison of the stagnation zone pressure ratio 3 p { 
in Tables 1 and 3 shows that the effect of the exit-plane choking is to provide (theoretically) up to 
8-percent higher pressures upstream of the model than for the Rayleigh calculation of the duct flow. 

To summarize the purely theoretical aspects of this work, the Rayleigh and Fanno calculations 
require similar auxiliary conditions for compatibility. The major difference between the two methods 
is that the Fanno duct-flow theory results in consistently higher pressure ratios throughout the entire 
range of up to 12.5 percent for all transmitted-wave solutions. 

Theoretical computations were also performed for an ideal isentropic model assuming a vena 
contracta with supersonic duct flow and no vena contracta with sonic duct flow. The computed 
transmitted-wave strengths are plotted in Figure 26 and discussed in a later section. 
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Experimental Results 


Qualitative Results 

The high-speed photographs provided by the framing camera for the majority of the test runs 
verified the existence of the wave phenomena predicted by the theoretical computations. For both 
models, at the incident-wave strengths of 4 to 5.5, the primary transmitted wave was accompanied by 
an unstable auxiliary wave which oscillated during the initial transient at a position about 6 to 7 hole 
diameters downstream of the model and then assumed the fixed stable position of a standing wave 
about 8 to 9 hole diameters downstream of the model. For incident-wave strengths of 8.5 and 
higher, the primary transmitted wave was accompanied by a left-facing traveling shock as was pre- 
dicted by the theoretical computations. In none of the test runs were any wave formations noted 
that could be interpreted as auxiliary rarefaction waves nor were there any pressure traces that showed 
the presence of rarefaction waves. 

Both models, in general, exhibited the same basic wave formations in the upstream reflected-wave 
regions and in downstream transmitted-wave regions with one exception. In the case of the model with 
0.078-in. holes, the primary transmitted wave rapidly developed as a plane wave and remained essen- 
tially a plane wave as it traveled downstream. The model that had the larger (0.250 in.) holes started 
as a curved wave with a parabolic bump at its center and very slowly flattened as it progressed down 
the tube. It did not lose its curvature completely while in the field of view and although the bump 
flattened, it was still evident at the edge of the optical field for the lower shock strengths. At higher 
incident-wave strengths, similar wave formations were noted with the appearance and disappearance of 
additional intermediate-wave formations between the primary wave and the auxiliary wave. The 
primary wave in traversing the length of the downstream window appeared to weaken and then 
strengthen and become a pair of plane traveling shocks, one located where the bump had been and one 
where the main curved front had been. The reason for the complexities associated with this model is 
a combination of effects that are present with the small hole model, but to a much lesser degree. 

The 0.078-in.-hole model had at least 68.5 percent of its frontal area involved in its core flow; 
therefore, as far as the optical system was concerned, the zone of Coanda flow at the channel walls 
was quite small. Recompression shocks at the wall affected only the edges of the primary wave which 
appeared to fully dominate the net flow as it progressively developed. In the case of the model having 
0.250-in.-diameter holes, only about 26.5 percent of the frontal area was involved in the central core 
flow, the remainder was associated with a Coanda flow at the channel boundaries. This permitted the 
core to progress more rapidly downstream producing the bump; the separated flow at the wall recom- 
pression produced the oblique shocks which then interacted with the core to provide the appearance 
of the curved wave and core bump. The intermediate waves that appeared and vanished as the primary 
wave traveled downstream could very readily have been caused by stagnation-region perturbations due 
to higher levels of turbulence in the field behind the reflected incident wave in the case of the larger 
holes. 

It is apparent that 0.250 in. was not an optimum choice of hole size because of this higher 
turbulence level and the domination of the field by wall effects. In addition, since only 17 holes were 
involved as compared to 1 79 for the first model (to provide the same axial porosity), there was less 
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likelihood of developing a uniform front as rapidly as would occur with the many refracting waves 
that interacted in the model with small holes. 

A model area scale factor of about 5 instead of the 10 used would have been preferable. This 
would have provided about 61 holes of 0.125-in. diameter to obtain the same porosity, but with much 
less turbulence and wall interference. 

Figure 20 shows several frames of the framing camera record from run 22. The incident-wave 
strength was 17.73 and the sequence shown is for that portion of the cycle after the incident wave 
reflected from the upstream face of the model. The reflected wave appeared to start out as a plane 
wave and rapidly bent at the edges when interacting with the boundary layer produced by the incident 
wave. A very turbulent field was produced behind the reflected wave. Downstream of the model, the 
primary transmitted wave was generated as a plane wave except at the region near the walls. The wall 
interaction produced the appearance of a thick blurred wave. It was rapidly followed by an auxiliary 
wave that was initially curved because of wall effects, but rapidly became plane at its central core 



Figure 20— Framing camera record for run 22 taken at 160 750 frames per 
second with incident-wave strength S/ = 17.73; arrows indicate true 
locations of model faces. 
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with classic lambda shocks produced by interaction with the boundary layer generated by the primary 
transmitted wave. Careful examination of the auxiliary wave indicated that the increased turbulence 
levels in the reservoir zone produced a progressive asymmetry of the wave. 

The increased turbulence was also evident in the oscillographs that were taken of the outputs of 
the pressure transducers located upstream and downstream of the model. Figure 21 shows the oscillo- 
graph records for run 5, which had an incident-wave shock strength of 5.27. The reflected-wave pres- 
sure oscillations were about twice that of the incident wave. The reflected wave had the typical abrupt 
increase in pressure initially, but then because of flow through the wall perforations and the blending 
of rarefactions from the holes and shocks from the solid portions of the wall, the pressure discontinuity 
became less abrupt as illustrated in Figure 21a. The downstream transmitted pressure waves appeared 
to have minimal turbulence. The records shown in Figure 22 were taken for run 22 at an incident- 
wave shock strength of 17.73. Increased turbulence levels were noted for the reflected waves. The 
pressure oscillations were of the order of five times those of the incident waves and the transmitted 
pressure waves were now very turbulent with the evidence of reinforcement by reflections. The reflected- 
wave heating effect on upstream transducer output was also apparent when the normal leveling of the 
pressure in trace was followed by a subsequent output increase at a lower rate. 

Figure 23 shows the similarity of traces for the large-hole model at an incident-wave shock 
strength of 21.86. The upstream turbulence was slightly higher than shown for the small-hole model 
in Figure 22, but the downstream turbulence was about doubled. In both cases the records became 
more difficult to analyze because of heating and pressure perturbations caused by the turbulence. 

Quantitative Results 

Tables 5 and 6 indicate the measured incident-wave speeds and incident-wave strengths for the 
primary and secondary models, respectively. 

Incident-wave speeds were measured over a 2-in. interval (2 in. from the model) from the oscillo- 
scope time-pressure histories and also over two 10-in. intervals upstream of the model by electronic 
counters. Plots of the three measurements as a function of position indicated that wave-speed 
attenuation measured by counter readings was random and, in general, negligible. The measurements 
taken from oscilloscope records provided the appearance of greater attenuation rates of incident speed 
than the counter measurements, but these differences were of the same order of magnitude as the 
probable error predicted by error analysis. Measurements from film records (e.g., Figure 20) were not 
usable for reasons discussed in the later section on error analysis. 

Therefore, the means of the counter readings were used to establish the shock -tube calibration 
curve shown in Figure 24. In general, most of the data points are within a ±5 percent error band. A 
small percentage of the plotted points exceed this limit, but are still within a 10 percent band. 

This is not unusual with this type of apparatus because, as has been indicated by Winkler, 1 varia- 
tions from ideal theory will occur because of variations in wave formation caused by anomalous 
diaphragm rupture, and partial mixing of the helium driver gas with the nitrogen driven gas. The 


1 E. Winkler, Naval Ordnance Laboratory, June 23, 1971. 
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Figure 21— Oscilloscope records for run 5; incident-wave strength Sj = 5.27 with primary model; calibration fre- 
quency = 10.071 kHz. (a) Upstream pressure traces taken at stations 1 18 and 120; calibration amplitude = 0.051 V/cm. 
(b) Downstream pressure traces taken at stations 2D and 4D; calibration amplitude = 0.050 V/cm. 



Figure 22— Oscilloscope records for run 22; incident-wave strength = 17.73 with primary model; calibration fre- 
quency = 9.955 kHz. (a) Upstream pressure traces taken at stations 118 and 120; calibration amplitude = 0.495 V/cm. 
(b) Downstream pressure traces taken at stations 2D and 4D; calibration amplitude = 0.049 V/cm. 



Figure 23— Oscilloscope records for run 40; incident-wave strength S f - 21.86 with secondary model; calibration fre- 
quency = 10.975 kHz. (a) Upstream pressure traces taken at stations 118 and 120; calibration amplitude = 0.502 V/cm. 
(b) Downstream pressure traces taken at stations 2D and 4D; calibration amplitude = 0.049 V/cm. 
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Table 5— Experimental values of incident-wave speeds and shock strengths for the primary model. 



Cathode-Ray Oscilloscope, 

Counters, 

Counters, 


Stations 1 1 8 and 1 20 

Stations 100 and 1 10 

Stations 1 1 0 and 1 20 

Run 








No. 

2^1 

M i 

U I 


u. 

Mj 

U i 




(fps) 


(fps) 


(fps) 

i 

4.23 

1.940 

2237 

— 

— 

2.030 

2340 

2 

5.35 

2.025 

2334 

— 

— 

2.220 

2553 

3 

— 

— 

— 

2.198 

2534 

2.213 

2552 

4 

5.31 

2.098 

2419 

2.228 

2568 

2.236 

2578 

5 

5.27 

2.120 

2441 

2.220 

2561 

2.240 

2581 

6 

5.53 

2.062 

2378 

2.239 

2581 

2.244 

2587 

7 

9.64 

2.810 

3240 

2.990 

3448 

2.995 

3454 

8 

9.73 

2.809 

3239 

2.994 

3452 

2.879 

3320 

9 

9.37 

2.664 

3072 

2.930 

3379 

2.929 

3378 

10 

8.50 

2.668 

3076 

2.785 

3211 

2.809 

3238 

11 

9.66 

2.895 

3338 

3.010 

3469 

3.030 

3496 

12 

11.91 

2.920 

3367 

3.224 

3717 

3.270 

3771 

13 

12.14 

3.170 

3660 

3.234 

3729 

3.266 

3766 

14 

10.71 

3.085 

3558 

3.077 

3548 

3.091 

3564 

15 

12.20 

3.083 

3555 

3.257 

3756 

3.291 

3795 

16 

11.96 

3.176 

3661 

3.259 

3757 

3.294 

3798 

17 

17.80 

3.567 

4113 

3.969 

4577 

3.965 

4571 

18 

17.78 

3.911 

4509 

4.155 

4791 

4.118 

4748 

19 

— 


— 

4.337 

5000 

4.054 

4674 

20 

17.08 

3.758 

4332 

4.079 

4704 

3.809 

4393 

21 

17.75 

3.896 

4492 

4.131 

4763 

3.816 

4400 

22 

17.73 

3.631 

4187 

4.041 

4660 

4.031 

4627 

23 

21.80 

4.423 

5100 

4.602 

5306 

4.531 

5225 

24 

20.99 

4.156 

4792 

4.553 

5249 

4.478 

5163 


greatest dispersion from ideal theory (not attributable to probable error) occurs at Mach numbers 
above 4. However, these are within the range reported by experimenters in well-established facilities 
(Reference 32). 

The data for runs where excessive leaks were known to exist were not included in Figure 24 as in- 
dicated in the later discussion of rest-region pressure errors. 

Values of reflected- and transmitted-wave strengths are tabulated as functions of the incident-wave 
strength in Tables 7 and 8 for the small-hole and large-hole models. All values represent the mean of 
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Table 6— Experimental values of incident-wave speeds and shock strengths for the secondary model. 



Cathode-Ray Oscilloscope, 

Counters, 

Counters, 

Run 

Stations 1 18 and 120 

Stations 1 00 and 1 1 0 

Stations 1 10 and 120 

No. 

2P1 

Mj 

(fps) 

M f 

(fps) 

Mj 

Uj 

(fps) 

25 

5.29 

2.063 

2378 

2.209 

2548 

2.223 

2563 

26 

5.16 

2.065 

2381 

2.192 

2527 

2.196 

2531 

27 

5.44 

2.200 

2537 

2.217 

2556 

2.254 

2599 

28 

5.07 

2.109 

2432 

2.161 

2492 

2.217 

2556 

29 

5.13 

2.021 

2330 

2.191 

2526 

2.203 

2540 

30 

9.61 

2.735 

3154 

2.952 

3404 

2.984 

3441 

31 

9.86 

2.663 

3071 

2.895 

3338 

2.923 

3370 

32 

9.57 

2.659 

3065 

2.921 

3368 

2.938 

3388 

33 

8.74 

2.797 

3227 

2.869 

3307 

2.875 

3315 

34 

9.23 

2.589 

2986 

2.801 

3230 

2.805 

3234 

35 

12.72 

2.815 

3245 

2.962 

3415 

2.995 

3454 

36 

1.215 

3.076 

3547 

3.072 

3542 

3.107 

3583 

37 

13.07 

2.889 

3331 

3.042 

3508 

3.098 

3572 

38 

13.52 

2.895 

3337 

3.038 

3503 

3.083 

3555 

39 

13.13 

2.980 

3438 

3.116 

3593 

3.119 

3600 

40 

21.86 

3.612 

4165 

3.591 

4140 

3.747 

4320 

41 

20.09 

3.760 

4336 

3.852 

4442 

3.820 

4405 

42 

21.50 

3.500 

4037 

3.690 

4255 

3.688 

4252 

43 

22.88 

3.500 

4036 

3.888 

4483 

3.676 

4239 

44 

— 

— 

— 

3.709 

4277 

3.705 

4271 

45 

19.48 

3.758 

4332 

3.824 

4409 

3.800 

4381 

46 

12.39 

2.747 

3167 

2.836 

3270 

2.855 

3291 


the peak values of two data traces exclusive of reducible heating effects. (Effects of pressure 
gage heating are discussed in a later section.) 

The reflected-wave strengths were plotted against incident-wave strengths for the two test 
models in Figure 25. Also plotted, for purposes of comparison, were the curves for the ideal 
solid-wall reflected-wave strengths and those values predicted by the Rayleigh and Fanno duct- 
flow theories. 

Since the curves for the theoretical data were less than 10 percent apart, the error bands 
shown are for +10 percent for the upper-limit solution and -10 percent for the lower-limit 
solution. For values of incident-wave strength below 15, the experimental data indicate that the 
simpler Fanno theory very closely predicts reflected wave strengths with an accuracy of better 
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Figure 24— Shock -tube calibration based on counter measurements of wave speeds. 


Table 7— Experimental values of mean peak reflected- and transmitted-wave strengths for the primary 

model. 


Run No. 

2 P 1 

3 P 1 

6P1 

Run No. 

2P1 

3P1 

6 P 1 

i 

4.23 

11.77 

2.35 

13 

12.14 

52.86 

7.49 

2 

5.35 

16.26 

2.88 

14 

10.71 

52.00 

6.59 

3 

— 

— 

— 

15 

12.20 

54.12 

7.33 

4 

5.31 

16.91 

2.66 

16 

11.96 

53.92 

7.39 

5 

5.27 

16.19 

2.97 

17 

17.08 

74.04 

9.77 

6 

5.53 

17.31 

3.11 

18 

17.78 

101.98 

10.53 

7 

9.64 

39.57 

5.66 

19 

— 

— 

— 

8 

9.73 

39.91 

5.73 

20 

17.08 

96.80 

10.62 

9 

9.37 

37.94 

5.42 

21 

17.75 

90.33 

8.91 

10 

8.50 

34.33 

5.01 

22 

17.73 

87.98 

10.05 

11 

9.66 

38.37 

5.26 

23 

2f.80 

115.74 

12.75 

12 

11.91 

55.27 

7.41 

24 

20.99 

126.10 

11.92 


Note: p 2 andp 3 measured at stations 118 and 120 ;p 6 measured at stations 2D and 4D. 
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Run No. 

2 P 1 

25 

5.29 

26 

5.16 

27 

5.44 

28 

5.07 

29 

5.13 

30 

9.61 

31 

9.86 

32 

9.57 

33 

8.74 

34 

9.23 

35 

12.72 


3 P 1 

iPl 

Rui 

16.23 

3.15 


14.34 

3.05 


16.48 

3.12 


15.52 

2.91 


15.81 

2.92 

i 

36.73 

5.26 

i 

39.23 

5.62 

i 

38.81 

5.48 

l 

35.96 

5.50 

i 

37.25 

5.67 

l 

55.90 

7.81 

i 



Table 8-Experimental values of mean peak reflected- and transmitted-wave strengths for the second- 
ary model. 


12.15 

13.07 

13.52 

13.13 

21.86 

20.09 

21.50 

22.88 


19.48 

12.39 


Note: p 2 and p 3 measured at stations 118 and 120;p 6 measured at stations 2D and 4D. 


IDEAL SOLID WALL (7 = 1-4) 

o RUNS WITH PRIMARY MODEL 
□ RUNS WITH SECONDARY MODEL _ 

FANNO DUCT FLOW 

RAYLEIGH DUCT FLOW 

+10% FANNO THEORY AND 

10% RAYLEIGH THEORY 

p = 5 psia T, = SSS 15 R 


55.14 

58.70 

57.90 

56.81 

110.17 

102.11 

106.35 

108.54 

97.65 

55.83 


7.71 

8.13 

8.20 

7.93 

11.67 

11.96 

12.85 

13.27 

13.34 

13.45 

8.11 


INCIDENT WAVE SHOCK STRENGTH S, = p 2 /p, 

Figure 25— Reflected -wave strengths as a function 
of incident-wave strengths. 
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than 5 percent. However, at values above an incident-wave strength of 15, the dispersion became 
significantly wide with several of the data points exceeding the +10 percent Fanno limit. This can be 
the result of the greater possibility of anomalous diaphragm rupture with the thicker diaphragm units 
required to produce the higher values of incident shock strength, the resultant mixing of driver and 
driven gas, and the region of very complex mixing of shock and expansion waves at the upstream face 
of the model by these gas mixtures. These complexities were indicated for very weak shock waves by 
Rudinger (Reference 13) for a single-hole pattern where incident pressure ratios of 1.0 to 2.4 were 
being considered. For the stronger waves and the multiplicity of perforations of concern here, it is 
reasonable to expect a wide range of deviations to be experienced. However, at no time did any of the 
data points fall outside the ideal solid-wall curve so that the ideal absolute limit on test data dispersion 
was not exceeded. 

The experimental values of reflected- and transmitted-wave speeds were tabulated for the small- 
hole and large-hole models in Tables 9 and 10, respectively, for the full range of incident-wave shock 
strengths. Although similar considerations concerning wave attenuation and correlation of measure- 
ments by oscilloscope and counter techniques apply here as for the incident-wave speed measurements, 
the close correlation of the reflected-wave speeds should be noted. This is not unusual when one 
considers that the reflected wave is slow moving and therefore can be more accurately measured from 
the oscilloscope record than either the incident or transmitted waves. Therefore, dispersion of values 
can be minimal. This is also seen to be the case for the lower values of the transmitted wave speeds. 
When incident shock strengths increase above a value of 10, the spread between oscilloscope and 
counter wave speeds also increases. 

It is interesting to note that the values of wave speeds shown in Tables 9 and 10 when compared 
with the values predicted by the Rayleigh and Fanno theories (Tables 2 and 4, respectively) indicate 
good correlation in terms of the predicted and experimental auxiliary-wave solutions as a function of 
the incident-wave shock strengths. 

Comparison of Theoretical and Experimental Results 

The average of the peak values of transmitted-wave pressure ratios as computed from measure- 
ments (at arrows in Figures 21 to 23) of the dual downstream transducer outputs were plotted as a 
function of the incident-wave strengths in Figure 26a. The theoretical curves for the auxiliary standing- 
and traveling-wave solutions of the Rayleigh and Fanno duct-flow models were plotted on the same 
graph. 

As predicted by theory and demonstrated physically by means of the high-speed photographic 
records, the data points for an incident-wave strength of about 5 fall within the predicted curves for 
the auxiliary standing wave. The data points for higher values of incident-wave strength coincide with 
the curves for the auxiliary left-facing traveling-wave solution. Figure 26b is an enlargement of 
the initial section of Figure 26a. It shows the major portion of the data points for both auxiliary 
solutions lie on the Rayleigh duct-flow curves. The theoretical curves are very close to each other for 
the same auxiliary-wave solution. The difference between theories for each solution is about 1 5 per- 
cent. The correlation of the experimental data with the Rayleigh solution, in general, is better than 
5 percent in the range of incident-wave strengths below 1 5 and the same order of correlation is 
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Table 9— Experimental values of reflected- and transmitted-wave speeds for the primary 

model. 




U R -u 

<2 ( f Ps) 


U rp (fpS) 


Run No. 

2 p l 

Cathode- 
Ray Oscilloscope, 
Stations 118, 120 

Counter, 

Stations 118, 120 

Cathode- 
Ray Oscilloscope, 
Stations 2D, 4D 

Counter, 
Stations 4D, 12D 

Counter, 

Stations 12D, 22D 

1 

4.23 

-929 

_ 

1845 

__ 

1801 

2 

5.35 

-1005 

-1008 

2032 

1953 

1972 

3 

- 

— 

-1018 

— 

1938 

1982 

4 

5.31 

-941 

-952 

2300 

2568 

2577 

5 

5.27 

-933 

-955 

1947 

*1976 

a 1976 

6 

5.53 

-925 

-940 

1930 

1969 

2025 

7 

9.64 

-1051 

- 

2650 

2612 

2640 

8 

9.73 

-1080 

- 

2325 

2616 

2644 

9 

9.37 

-1062 

- 

2581 

2560 

2585 

10 

8.50 

-1044 

- 

2396 

2444 

2473 

11 

9.66 

-1102 

- 

2623 

2631 

2658 

12 

11.91 

-1203 

- 

2907 

2907 

2946 

13 

12.41 

-1171 

-1169 

2693 

2899 

2934 

14 

10.71 

-1087 

-1070 

2508 

2757 

2793 

15 

12.20 

-1117 

-1132 

2911 

2925 

2957 

16 

11.96 

-1196 

-1189 

2562 

2925 

2952 

17 

17.80 

-1413 

-1417 

3224 

3440 

3474 

18 

17.78 

-1447 

-1465 

3487 

3592 

3631 

19 

- 

- 

-1432 

- 

- 

- 

20 

17.08 

-1409 

-1385 

3470 

a 3605 

a 3605 

21 

17.75 

-1424 

-1394 

3472 

3552 

3623 

22 

17.73 

-1413 

-1422 

3351 

3480 

3544 

23 

21.8 

-1552 

-1572 

3707 

a 3921 

a 3921 

24 

20.99 

-1559 

-1545 

3670 

3872 

3913 


a These values axe computed from readings at stations 4D and 22D. 


evident at higher shock strengths with a few isolated points just outside the ±10 percent error band of 
the Rayleigh theoretical curves. Those values outside the error band were dispersed equally about the 
theoretical curve and can be attributed to the greater influence of leakage on the computation of 
transmitted-wave strengths and the influence of the complex transient processes as a result of the more 
turbulent mixing of shock and expansion waves during downstream transients. Since the data points 
involved were those obtained with the secondary model (0.250-in.-diameter holes), this greater devia- 
tion can also be attributed to the very complex three-dimensional fields produced by a combination of 
the normal downstream transient fields coupled with the significant Coanda effect of the wall inter- 
actions. These very turbulent and complicated flows can very readily deviate from the simple one- 
dimensional theories utilized because of greater variations in entropy and localized shock phenomena. 

Error Analysis of Experimental Data Systems 

Pressure Instrumentation 

Each of the components of the pressure data system affects the accuracy of the acquired time- 
pressure histories. To realistically assess overall reduced data limitations, the transducer and its 
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Table 1 0— Experimental values of reflected- and transmitted-wave speeds for the secondary 

model. 




Ur -« 2 < f P s > 


U T { fps) 




Cathode- 

Counter, 

Stations 118, 120 

Cathode- 



Run No. 

2Pl 

Ray Oscilloscope, 
Stations 118, 120 

Ray Oscilloscope, 
Stations 2D, 4D 

Counter, 
Stations 4D, 12D 

Counter, 

Stations 12D, 22D 

25 

5.29 

-925 

_ 

2074 




26 

5.16 

-949 

-953 

2032 

1998 

1979 

27 

5.44 

-912 

-983 

2136 

2043 

2035 

28 

5.07 

-919 

-925 

1959 

1974 

2025 

29 

5.13 

-925 

-929 

2037 

2011 

2006 

30 

9.61 

-1024 

- 

2711 

2619 

2623 

31 

9.86 

-1006 

- 

2914 

2572 

2571 

32 

9.57 

-1069 

- 

2638 

2591 

2600 

33 

8.74 

-1028 

- 

2648 

2544 

2553 

34 

9.23 

-1059 

- 

2524 

2486 

2491 

35 

12.72 

-1053 

-1062 

2682 

a 2717 

a 2717 

36 

12.15 

-1148 

- 

2814 

2776 

2835 

37 

13.07 

-1070 

-1072 

2747 

2774 

2910 

38 

13.52 

-1092 

-1087 

2749 

2773 

2823 

39 

13.13 

-1135 

-1137 

2750 

2770 

2813 

40 

! 21.86 

-1090 

-1097 

3796 

4201 

4861 

41 

20.09 

-1330 

-1331 

3436 

3389 

3418 

42 

21.50 

-1286 

-1301 

3148 

a 3688 

a 3688 

43 

22.88 

-1272 

- 

3060 

3335 

3556 

44 

- 

- 

-1411 

- 

3381 

3511 

45 

19.48 

-1376 

-1459 

3082 

3389 

3443 

46 

12.39 

-1010 

-1133 

2304 

2473 

2670 


a These values are computed from readings at stations 4D and 22D. 


auxiliary electronics, the recording instrumentation, the calibration trace instrumentation, and data 
reduction methods must be considered in terms of their influence on pressure amplitude and time base 
errors. This is especially important since wave speeds deduced from pressure traces should correlate 
with those measured directly by electronic counters. 

The amplitude errors due to combined calibration errors of standards, transducer, and charge 
amplifier are no greater than 1.2 percent of data point value. The error due to calibration trace ampli- 
tude instability can be as high as 4 percent at 50 mV to 1 percent at 500 mV. The data reduction 
errors based on a 0.005-in. trace resolution produces an additional error of 5 percent for traces of 
0.10-in. length to 1 percent for traces of 0.500-in. length. This means that for low-level pressure sig- 
nals, it is conceivable that cumulative pressure magnitude errors of 10 percent maximum are possible. 
This error is reduced for the higher values of pressure to approximately 3.5 percent for reduced data. 

Since the pressure traces were also used to obtain incident- , reflected- , and transmitted-wave 
speeds, the stability and linearity of the oscilloscope sweep rate, the calibration trace frequency, and 
the trace reading error become significant factors. The oscilloscope manufacturer’s data (References 
33 and 34) indicate a maximum full-scale sweep error of 3 percent at any point along the central 8 cm 
of the 10-cm graticule. At the frequencies set for the pulse generator that provides the calibration 
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square wave, a !4 percent rate instability was recorded during the test program. Therefore, equipment 
errors alone can provide a maximum error of 314 percent. The trace reading error for low-velocity 
waves (such as reflected waves) results in reduction errors of about 1 percent; however, for fast-moving 
waves where the trace length between points of interest are on the order of 0. 100 in., the reduction 
error can increase to a maximum of 5 percent. Therefore, if the pressure traces are used to establish 
local wave speeds, the possible error range can be from 414 to 814 percent of the true value. 

The final factor that must be considered is the effect of transducer heating errors due to the high 
stagnation temperatures occurring behind the reflected incident shock at the upper end of the incident- 
wave test range. 

Temperatures on the upstream side of the model can increase by 3000 F° over the initial ambient 
of 535° R. Since the transducer manufacturer claims a 0.01 percent output for every Fahrenheit- 
degree temperature rise, if the transducer assumed the gas temperature, it could conceivably have an 
error of + 30 percent. However, the gas medium is transparent to radiation in the temperature range 
expected, and the testing time is too short for the transducer to establish these temperatures through 
conduction and convection heat transfer. 

Discussions of this problem with personnel of the Naval Ordnance Laboratory 2 indicated that they 
have considered this problem on a limited basis. Aronson of the Aeroballistics Division has attempted 
to insulate such transducers by means of 6-mil black electrical tape. He has measured amplitude 
differences of 3 to 5 percent for incident-wave Mach numbers of 5 to 8 by this method. The effect of 
additional thicknesses on reducing thermal error were not considered because of a falloff in dynamic 
response with the 6-mil tape. This technique was established on a trial and error basis because of the 
complexity of an analytical approach to the heat transfer involved. Therefore, there is no way to 
establish the absolute effectiveness of this technique. McLain (Reference 35) has also noted piezo- 
electric transducer baseline shift during combustion stability studies of N 2 H 4 /N 2 0 4 reactors. He 
attributed this shift to radiation from the combustion products and assessed the amplitude of the shift 
by a composite time-pressure history lasting 5 ms with a 100 jus/cm sweep rate. 

The effect of heating was taken into account in the current investigation by careful examination 
of the pressure trace. Unlike McLain’s situation, the time-pressure traces were relatively clean and it 
was possible to visually determine where the maximum output occurred prior to heating by the level- 
ing off of the trace and then by the abrupt increase in slope with additional amplitude growth as heat- 
ing occurred. This effect is sufficiently distinctive so that it is not obscured by the high turbulence 
levels occurring after reflection from the perforated wall of the model. In Figure 22a, this effect is 
plainly evident in the upper trace of the reflected wave. The peak indication is marked by a white 
arrow; however, this includes the heating effect and is followed by the cooldown portion of the trace. 
Exactly 3 cm to the left of the white arrow (vertical lines on graticule are 1 cm apart) is the place 
where peak pressure was actually measured. It is seen that at this location in time there is a zero slope 
to the mean pressure indicating that the trace had hit the limit of its initial rapid rise. Further in- 
crease in amplitude was attributed to the slower rate occasioned by heat transfer. The lower trace on 
this same photo reaches its peak pressure 2 cm to the left of its corresponding white arrow. All traces 
produced by incident-wave strengths above 1 5 were reduced in this manner. The traces at lower 
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values of wave strength were not sufficiently distinctive to employ this visual technique, and therefore 
temperature was not considered to be of consequence in contributing to data reduction errors. 

Rest-Region Pressure Errors 

To very accurately reproduce initial conditions, a master vacuum gage having 1.0-mm-Hg divisions 
was used to charge the driven section to 501.5 mm Hg with nitrogen. Unfortunately, the gage had to 
be valved off prior to firing the shock tube to prevent damage to it by reflected shocks transmitted in 
the lines. Therefore, for every run there was a period of about 2 min before a firing in which the rest- 
region pressure was not being monitored with precision equipment. 

During the testing of the secondary model, leaks developed in the driven section after run 34 so 
that it was not possible to pump down on the scavenging part of the test cycle to the required 0.5 mm 
Hg. Attempts were made to compensate for the leak by a vacuum bleed during test countdown. 
Examination of the data for runs 35 through 46 showed that shock strengths computed on the basis 
of the rest-region pressure being p x = 5 psia resulted in incident wave-speed values considerably lower 
than predicted by ideal theory. On the basis of ideal theory, it was computed that the rest-region 
pressures ranged from an average value of 6.01 for runs 35 to 39 to 6.66 for runs 40 to 46. Apparently 
the leak had grown progressively worse with further testing. To avoid complete loss of the data, 
those presentations that involved the pressure ratio in only one of the variables excluded the question- 
able data points; those presentations for which the error was compensated by having the rest-region 
pressure common to the dependent and independent ordinate retained all data for purposes of data 
analysis. 

Wave-Speed Instrumentation Errors 

Although the electronic counter system accuracy depends primarily on the precise adjustment of 
its crystal clock in accordance with the U.S. Bureau of Standards worldwide radio time signal (WWV) 
standards, the auxiliary circuitry used for the amplification and the gating of timing signals could 
introduce considerable error. Since differences between individual readings were utilized, the zero 
errors for common stop signals were of importance. These produced no difficulty on forward travel- 
ing waves; but on reflected-wave circuits, zero errors ranging from 5 to - 10 ms were noted. These 
errors were a function of wave strength, wave shape, and the aging characteristics of the equipment. 

For incident-wave strengths in the range of 9 to 10, it was not possible to obtain reliable triggering of 
the pulse-counting gate so that reflected-wave speeds were not obtained directly by counter measure- 
ment. On other runs, where required, correction factors to compensate for zero error were included. 

Optical System Errors 

The optical instrumentation, although adequate for qualitative studies, provided data of ques- 
tionable validity when considered quantitatively. There seemed to be excessive, random deviations 
from the mean of wave velocities measured from the framing camera film strips. The measuring errors 
at low velocities were quite high since the frame-to-frame wave movement was on the order of 3 to 
10 times the smallest measurable increment of 0.010 in. Thus errors on the order of 10 to 30 percent 
were possible if the optical components were perfect and introduced no additional inaccuracies. At 


48 



higher wave speeds, the error was on the order of 5 to 10 percent; the maximum and minimum devia- 
tions from the mean were significantly less for these runs, but still not comparable to data obtained 
from counter or oscilloscope records. 

This was evident when velocity-time measurements obtained from the film frames were plotted. 
Very random oscillatory values for wave speed were obtained. The oscillations were of a significantly 
high percentage of the velocity amplitude. In addition to the errors in resolution were the distortion 
errors resulting from the plastic windows that were substituted for damaged optical flats. They were 
out of flat and parallel by as much as 0.005 to 0.010 in. and had considerable waviness and local 
discontinuities (such as rounded edges, pits, fine scratches) that served to obscure weak wave fronts 
and produced interference bands at model/window interfaces. 

The camera optics also contributed to deviations from an ideal system in the space dimension. 
Time resolution was more than adequate since mirror turbine speeds were continuously monitored by 
an electronic counter built into the camera control console. The major error sources were the pre- 
cision of the alignment of the individual relay lens and mirror systems (80 of them) and the variable 
length of the optical lever arm from frame to frame. It is not possible to attach numbers to these 
errors since the alignment was accomplished visually on each system. A rigorous study of the vibration 
modes of the 80 mirror mounts, the relay lens systems, and the input optics would probably indicate 
the possibility of deviations due to mechanical and aerodynamic excitation of elastic elements of the 
component structures and possible magnification of these vibrations due to transmissibility of the 
camera housing and base supports. 

Because of the large number of intangibles that did not permit an adequate quantitative error 
analysis, the wave speeds measured from the camera data were not presented. However, except in the 
immediate vicinity of the model (about 1/8-in. upstream and downstream), the high-speed record 
helped resolve and verify the adequacy of the assumption of one dimensionality of the wave fronts, 
the verification of auxiliary transmitted-wave solutions, and extent of wall effect on channel flows. 

Test Medium Variability Effects 

The test fluid was dry nitrogen with a dewpoint of - 105° F; for all practical purposes the presence 
of trace quantities of water vapor was negligible. However, during runs 35 through 46, ambient air 
leaked into the rest region. Since the ambient atmosphere had a 50 percent relative humidity, compu- 
tations were performed to examine the effect of air and water vapor contamination on the nitrogen 
charge of the driven section. The effect on y was negligible; for final rest-region pressures of 6 and 6.5 
psia in the driven section, the gas constant was reduced by 0.24 and 0.34 percent, respectively. Since 
the major influence of the gas constant is reflected in the computation of the sound speed, which is a 
function of the square root of the gas constant, the error contributed by gas constant change can also 
be neglected. 

CONCLUSIONS 

Limitations of Theoretical Work 

The theory as outlined herein predicted the reflected-incident-wave strengths to within a 5 percent 
accuracy for incident-wave strengths of 4 to 1 5; for incident-wave strengths between 15 and 25, the 
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reflected-wave strengths were within 10 percent accuracy. These values are valid when the Fanno duct- 
flow model is employed. 

When using the Rayleigh model for predicting reflected-wave strengths, the predicted values were 
on the low side by slightly over 5 percent for the full range of incident shock strength. 

The Rayleigh theory verified the nature of the transmitted auxiliary waves for the full range of 
the incident shock strengths considered. It also predicted transmitted-wave shock strengths with an 
accuracy better than 5 percent in the range of incident-wave strengths of 4.5 to 15, and predicted 
transmitted-wave shock strengths with an accuracy of ±5 to +10 percent in the incident-wave strength 
range of 1 5 to 25. 

The Fanno theory predicted the nature of the transmitted auxiliary waves more accurately than 
the Rayleigh theory; but this simpler approach enabled the prediction of transmitted-wave strengths 
to no better than 10 to 15 percent over the entire range of incident-wave shock strengths. 

Examination of the data in Figure 26 indicates that the transmitted-wave shock strengths were 
60 percent of the incident-wave strengths throughout the test range. Transmitted-wave speeds averaged 
about 77.5 percent of incident-wave speeds. 

Limitations of Experimental Work 

The experimental program was hampered by a less than optimum optical system. The use of 
plastic windows degraded optical performance to the degree that quantitative data were not reliable. 

The defects (previously outlined) also served to obscure details of the flow in close proximity to the 
model surfaces. However, these same windows were adequate for the purpose of qualitatively defining 
physical details of wave formulation, boundary-layer interaction, and wall effects. 

Piezoelectric transducers performed satisfactorily in providing time-pressure histories with rest- 
region conditions as the zero ground level. In perfectly sealed systems, a precision vacuum gage read- 
ing of rest-region conditions provided data of desirable accuracy. However, leaks can occur during 
test runs. Therefore, the vacuum gage should be supplemented by a rugged transducer or suitable 
settling chamber to enable continuous monitoring of the rest-region conditions. 

For wave speeds of about 1500 ft/s, pressure gage traces were not a sufficiently accurate means 
of assessing wave speeds for the 2-in. gaging distance utilized. If transducers are spread to a 6-in. 
spacing, then the probable timing accuracy can be guaranteed at least to ±3 percent of the oscilloscope 
sweep period. 

Counter readings of wave speeds were preferable in all cases except for reflected waves, where the 
pulse-counting gates apparently were sensitive to spurious triggering by the very turbulent fields and 
resulting unstable wave forms that accompanied the reflections. For measurements of these wave 
speeds, the analysis of the paired transducer outputs was the most reliable method. 

Recommendations for Additional Work 

The following recommendations are suggested for the purpose of supplementing and amplifying 
the theoretical and experimental work described herein. 
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Additional theoretical work would be to examine the effect on qualitative and quantitative cor- 
relations of the various duct-flow theories with the following variations: 

(1 ) Assume a constant y and Z in the Rayleigh method. 

(2) Allow for boundary-layer displacement thickness in computing expansion ratio at the duct 
exit plane in the Rayleigh method. 

(3) Recompute the Fanno duct-flow theory for a variable y and Z. 

(4) Average the optimum combination of Rayleigh and Fanno methods to provide the initial 
conditions for an influence coefficients solution of the problem. 

(5) Evaluate all methods indicated above as to time required and relative quality of results. 

(6) Consider a broader range of porosity, duct to channel area ratio, and L/d. 

Additional experimental work would include significant improvements in experimental techniques 
to enable the reduction of the data system’s contribution to experimental uncertainty. The following 
suggestions are offered to indicate considerations that could prove fruitful in accomplishing this goal: 

( 1 ) Replace the plastic windows with glass optical flats. 

(2) Consider an instrumented technique of focusing individual mirror systems in the framing 
camera. 

(3) Consider the addition of rigid sectors for mounting adjustment screws to lock mirror 
positions and eliminate vibration of mirror mounts. 

(4) Rigidize the camera support structure. 

(5) Investigate methods of reducing pressure transducer heating while not degrading response 
and also methods of evaluating the effects. (Consider the use of thin-film ceramic and 
reflective coatings, etc.) 

(6) Investigate ways of keeping the rest-region precision pressure gage in the system for static 
measurements that would isolate it from shock pressures when the tube is fired. 

(7) Consider the use of perforated models to study the dynamic growth of boundary layers in 
shock tubes by the interaction of the auxiliary transmitted waves with the growing layer 
left behind by the primary wave. 

(8) Conduct pressure-time-position studies to more critically evaluate wave field development. 

(9) Include thin-film gage temperature instrumentation to obtain temperature-time-position data. 

It is conceivable that the consideration of methods and implementation of these methods to 
accomplish these recommendations would enable improved correlation between theory and experi- 
ment in future work and would indicate other fruitful areas for fundamental investigation. 


Goddard Space Flight Center 

National Aeronautics and Space Administration 
Greenbelt, Maryland December 14, 1971 
160-20-51-12-51 
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